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ABSTRACT 

The  reduction  of  structural  vibrations  in  ships  and  submarines  is  a  long  standing 
concern  of  the  Navy.  Waveguide  absorbers  are  very  effective  devices  which  can  be  ap- 
plied to  this  problem.  This  study  evaluates  the  increase  in  vibration  damping  of  a  plate 
structure  across  a  broad  frequency  range  using  light  weight  beam  waveguide  absorbers. 
Viscoelastic  and  constrained  layer  beam  waveguide  absorbers  were  studied  both  theore- 
tically and  experimentally.  Impedances  of  the  waveguide  absorbers  at  the  attachment 
point  were  predicted  using  both  Bernoulli-Euler  and  Timoshenko  beam  theory  for  the 
viscoelastic  beam  and  using  sixth  order  beam  theory  for  the  constrained  layer  beam.  The 
theoretically  predicted  impedances  were  compared  with  the  experimental  measurements. 
Results  from  random  vibration  tests  of  a  plate  structure  showed  significant  increases  in 
damping  over  a  broad  frequency  range  ( 100  Hz  -  2  KHz)  when  the  waveguide  absorbers 
were  attached  on  the  plate. 
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I.     INTRODUCTION 

Suppressing  noise  and  vibrations  of  ships  and  submarines  is  very  important  for  the 
Navy  and  has  been  one  of  the  Navy's  long-standing  concerns.  Much  work  has  been  done 
on  the  development  of  vibration  control  means  such  as  isolation,  detuning,  viscoelastic 
damping  and  dynamic  absorption.  All  of  these  approaches  have  been  studied  rather 
thoroughly  including  their  range  of  applicability  and  definite  limits  on  their  vibration 
and  noise  reduction  capability  within  acceptable  weight  and  volume  increases. 

Recently,  a  waveguide  absorber  concept  has  been  developed  [Ref.  lj.  which  can 
provide  simple  and  highly  effective  vibration  control  over  a  wide  frequency  range.  A 
"Waveguide''  is  a  structure  along  which  vibrational  waves  can  travel.  If  one  end  of  a 
waveguide  is  attached  to  a  vibrating  structure,  some  vibration  energy  will  travel  along 
the  waveguide  in  wave  forms.  If  the  waveguide  is  treated  with  a  high  energy  dissipation 
scheme,  the  damping  of  the  waveguide  causes  the  amplitude  of  the  waves  to  decrease 
as  they  travel  and  the  waveguide  may  be  expected  to  remove  vibrational  energy  from  the 
structure. 

In  previous  studies  by  Ungar  and  Kurzweil  [Ref.  2]  and  by  Ungar  and  Williams  [Ref. 
3],  semi-infinite  beams  and  exponentially  tapered  semi-infinite  beam  waveguide  absor- 
bers were  studied.  However,  theoretical  prediction  of  the  driving  point  impedances 
showed  wide  discrepancy  from  the  experimental  results  since  the  experiments  were  per- 
formed with  finite  length  beams  and  the  effects  of  viscoelastic  materials  on  beams  were 
not  considered  in  theoretical  calculations. 

In  this  investigation,  two  kinds  of  high  damping  beam  waveguide  absorber  are 
subjected  to  study:  viscoelastic  beam  and  constrained  layer  beam  waveguide  absorbers. 
The  driving  point  impedances  are  theoretically  predicted  for  the  finite  length  beams  of 
highly  viscoelastic  behavior.  Experimental  studies  are  also  performed  on  the  impedances 
of  these  waveguide  specimens  and  on  their  contribution  to  damping  increase  of  a  test 
plate  structure. 


II.     THEORETICAL  ANALYSIS 

A.     LOSS  FACTOR  CONTRIBUTION  OF  A  WAVEGUIDE  ABSORBER 

Consider  a  linear  structure.  5,  under  a  harmonic  load 

F=fo  sin  ujt  (2.1) 

at  location  c  (Figure  1  on  page  3).  The  structure  5  will  vibrate  with  the  energy  of  vi- 
bration 

W^\MV2m  (2.2) 

where  M  represents  the  total  mass  of  the  structure  and  Vm  is  the  magnitude  of  average 
velocity.  The  energy  dissipated  in  5  per  cycle.  Da  depends  on  the  original  loss  factor.  >/„. 
of  the  structure  and  the  following  relation  holds: 

>/o  =  -Trrr-  (2-3) 

For  a  lightly  damped  structure.  vj0  is  much  smaller  than  1  and  D0  is  much  smaller  than 
2nlV0.  If  a  point  a  of  a  second  body  B.  a  waveguide  absorber,  is  attached  to  body  S 
at  location  b  of  body  S  .  internal  force,  Fb,  will  interact  between  body  S  and  body  B  and 
some  vibration  energy,  D  will  be  dissipated  in  body  D  (Figure  1  on  page  3).  Moreover, 
this  internal  force  influences  the  vibrational  motion  of  body  S  and  may  reduce  the  en- 
ergy of  vibration.  IV.  of  body  S  under  the  same  harmonic  load  at  the  same  location  c  . 

1  he  interaction  between  bodies  5  and  D  depends  on  the  impedances  of  bodies  .S  and 
D  at  the  attachment  point,  Z,  and  Zb.  Impedance  is  the  ratio  of  the  harmonic  force,  F. 
acting  on  a  point  a  of  the  structure  to  the  velocity.  V,  of  the  point. 

Z--£f  (2.4) 

and  can  be  describd  by  a  complex  number, 

Z=R  +  iX  (2.5) 
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Figure  1.      Vibrating  structure  with  waveguide  absorber. 


where  R  is  the  real  part  of  Z  and  resistance  component,  or  simply  resistance,  and  X  is 
the  imaginary  part  of  Z  and  reactive  component,  or  reactance. 

The  energy  loss  per  cycle,  D,  of  a  structure  by  an  attached  waveguide  absorber  is, 
as  shown  by  the  previous  study  [Ref.    2  ,  and  Ref.  3], 
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where  co  denotes  the  vibration  frequency  in  radian.'sec,  Vs  represents  the  velocity  of  the 
structure  at  the  attachment  point  of  the  structure  before  the  attachment  of  the  wave- 
guide absorbers.  Zh  and  Zr  denote  the  impedance  of  the  waveguide  absorber  and  the 
structure  at  the  attachment  point,  and  Rb  represents  the  real  part  of  the  impedance  Zb. 

From   the   equations    2.2,    2.3    and    2.6,    the   energy   loss   factor.   ;/,    due    to    the 
contribution  of  an  attached  waveguide  absorber  can  be  expressed  as  follows: 


R,       '   Vm  ' 
'7=7777 7 l2-7^ 

J 

Equation  (2.6)  shows  that  the  vibration  energy  loss  of  a  structure  through  the  wav- 
eguide absorber  depends  on  the  magnitude  of  the  velocity,  Vr,  at  the  attachment  point 

on  the  structure.  If  this  attachment  point  is  a  node  for  which  V,  =  0.  then  the  enersv  loss 

V. 
{actor.  ;/.  equals  zero.  II  the  attachment  point  is  an  antinode.  then  |  -p-  |2  can  be  rela- 
tively large.  Equation  (2.7)  indicates  that  the  loss  factor  contribution  is  small  if  the  im- 
pedance of  the  beam.  Zb  ,  is  considerably  greater  than  the  impedance  of  the  structure. 
7...  If  Z,  is  much  smaller  than  Zs.  then  the  structure  will  not  be  affected  by  the  absorber 
and  the  vibration  energy  of  the  structure.  \\\  in  equation  (2.2),  will  not  be  reduced  much. 
Therefore,  impedance  matching  between  the  structure  and  the  waveguide  absorber  is 
important  to  get  the  highest  energy  loss  factor, ;/.  So.  analytic  methods  which  can  predict 
the  impedance  of  beams  using  the  theory  of  elastic  wave  propagation  in  beams  are  de- 
veloped in  the  next  section. 

B.     IMPEDANCE  OF  BEAM  WAVEGUIDE  ABSORBERS  FROM  WAVE 
PROPAGATION  THEORY 

1.     VISCOELASTIC  BEAM  USING  BERNOULLI-EULER  BEAM  THEORY 

Previous  research  [Ref.  3]  considers  the  impedances  of  the  infinite  viscoelastic 
beam  using  Bernoulli-Euler  beam  theory.  In  this  section,  this  previous  work  is  expanded 
to  include  the  finite  viscoelastic  beam.  The  viscoelastic  beam  will  be  excited  harmonically 
at  its  center  by  a  transverse  force,  F  as  shown  in  Figure  2  on  page  6.  For  a  beam  which 
is  free  of  lateral  loading  and  under  an  assumption  that  cross-sectional  areas  remain  plane 
and  normal  to  the  neutral  axis,  the  equation  of  motion  becames  [Ref.  2]. 

a4  a2 

EKI—^-  +  pA-^r  =  Q  (2.8) 

OX  ct 

where  yipcj)  is  the  transverse  displacement,  and  Ex  is  the  complex  moduli  of  the  beam: 

£x  =  £(l  +  /^ve),  (2.9) 

and. 

E  :  modulus  of  elasticity  of  the  beam 


rjve  :  loss  factor  of  the  viscoelastic  material 
m  :  mass  per  unit  length  of  the  beam 
p  :  density  of  the  beam  material 
A  :  cross  sectional  area  of  the  beam 
I :  moment  of  inertia  of  the  cross  sectional  area 

The  general  form  of  >(.w)  for  a  propagating  harmonic  wave  in  this  Bernoulli- 
Euler  beam  is 

y(x,t)=  Y0eKkXx-Mt)  (2.10) 

where  kx  is  the  complex  wave  number,  and  co  is  the  frequency  of  the  propagating  wave. 
Inserting  the  wave  equation  (2.10)  into  the  equation  of  motion  (equation  2.S).  the  dis- 
persion relation  is  obtained  as  follow: 

ExI(kxf  =  p.ho2  (2.11) 

This  equation  gives  the  four  different  complex  wave  numbers  as  function  of  frequencies, 

-l 


*r=  ik(jl  +  vVe)    '\cos( -r-^-ZskK t-^-)J  (2A2.a 


4 


v-i/4r       ,  tan   lVve  ,      .   •   ,  tan   \]ve 


/<2X  =  -/7c(v  1  +  n\e )  xl*[cos( J**- )  -  /  sin< -^- )]  (2. 1 l.b) 


v-i/4,      t  tan  V«  s      .   .      tan   '?/v. 


A;  =  -  k(^'  1  +  ^  )-*'>«( ~^)  -  /  sin( j-11- )]  (2.12.C) 


,-1/4,       ,  tan   ''/v.?  ,  ,   t™   ';/, 


£4X  =    A-(v  1  +  n;,  r'\cos{ -p^- )  -  i  sm( -^L  )]  (2. 1 2.,/) 

where  subscript  ve  is  an  abbreviation  of  Viscoelastic  beam  using  Bernoulli-Euler  beam 
theorv  and  k  is  the  wave  number,  given  bv 

4  _    (Q_m    _    CO_£A    _  f    CO      i  n  ,  ,s 

a-  denotes  the  radius  of  gyration  of  the  cross  sectional  area,  (x/~r)  and  c,  denotes  the 
longitudinal  wavespeed  in  the  beam  material  {\j ~t)-  Therefore,  the  equation  (2.10)  is 
expanded  using  wave  number  k\,  /<2\  k;  and  k$> 


Figure  2.       Coordinates  and  sign  conventions  for  a  semi-infinite  (finite)  viscoelastic 
beam. 


y(xj)  =  (i^e   '    +  r2<?   -    -I-  r3<?  3    +  >at?       )e 


(2.14) 


where  ?",,  }\,  Y3  and  F4  are  constants  whose  values  depend  on  the  boundary  conditions 
and  x  is  the  axial  coordinate.  The  complex  transverse  force,  F,  is  denoted  by 


F(x)  =  F7 


67'-V(.v) 
c/.V3 


(2-15) 


a.     /AT/AY  FE  VIS  CO  EL  AS  TIC  BE.  IM 

Let's  consider  the  infinite  viscoelastic  beam.  IfyCr)  is  to  remain  finite  for  a 
large  x,  then  }\  must  vanish.  If  there  is  no  wave  coming  toward  the  origin  from  the  po- 
sitive x  direction,  then  Y2  must  vanish.  Yx  and  }\  in  equation  (2.1-4)  are  determined  from 
two  boundary  conditions  at  x  =  0,  as  functions  of  j/  and  6,  as  follows: 


y(Q)=y 


(2.16) 


&&-B  (2.17) 

ax 

At  the  driving  point,  the  input  transverse  force,  F(0),  can  be  calculated  using 
the  equation  (2.15),  which  can  be  expressed  using  impedance  and  velocity  at  x  =  0  as 
follows: 

F(0)  =  Zx  V{0)  +  Z2Q(0)  (2.18) 

The  complex  transverse  velocity,  F,  and  the  complex  angular  velocity,  LI,  are  related  to 
the  corresponding  displacement,^  and  6  of  equations  (2.16  and  2.17).  as  follows: 

V       P 

i-  =  ir  =  -^  (2.19) 

Z,  represents  the  impedances  of  the  transverse  force  at  the  driving  point  of  the  center 
of  the  infinite  viscoelastic  beam.  If  there  is  no  rotation  at  the  center  of  the  beam,  that 
is  Q(p)  =  0.  the  equation  (2.18)  becomes 

F(0)  =  ZxV(0)  (2.20) 

Therefore,  the  impedance  of  an  infinite  Bernoulli-Euler  beam  at  the  driving  point.  Zivt, 
can  be  calculated  from 

Z^  =  -^~T  (2-2D 

where  subscript  ive  is  an  abbreviation  of  the  Infinite  Viscoelastic  beam  using  Bernoul- 
li-Euler beam  theory. 

b.     FINI TE  VISCOELA  S  TIC  BE  A  A I 

For  a  finite  beam,  four  constants  Y\,  Y2,  Y3  and  K4,  are  determined  from  four 
boundary  conditions  at  x  =  0  and  x  —  /,  as  follows: 

y{0)=y  (2.22) 

dy{0) 


dx 


=  6 


">  ■>" 


.? 


■£x/— ^-V-  =  ()  (2.24) 

</x2 


,  A(/) 

El— '—-  =  <)  (2.25 

dx3 

From  the  same  considerations  as  those  in  the  previous  section,  the  impe- 
dance ZA„ 

Ffvp{0)    i 

is  determined  as  function  of  frequencies,  where  the  subscript  fve  is  the  abbreviation  of 
the  Finite  Viscoelastic  beam  using  Bernoulli-Euler  beam  theory. 

2.     FINITE  VISCOELASTIC  BEAM  USIiNG  TIMOSHEiNKO  BEAM  THEORY 

The  Bernoulli-Euler  beam  theory  was  obtained  using  two  assumptions  which 
neglect  shear  deformation  and  rotational  inertia.  The  Bernoulli-Euler  beam  theory  may 
be  useful  for  low  frequencies  and  long  wavelengths.  However,  the  Timoshenko  beam 
theory  gives  more  accurate  simulation  at  wider  frequency  ranges  and  wave  lengths. 
Therefore  we  will  analyze  the  impedances  of  the  finite  viscoelastic  beam  using  Timosh- 
enko beam  theory  in  this  section.  Two  independent  variables,  the  transverse  displace- 
ment, y{x,t),  and  angular  displacement.  (j){x.i).  are  defined  in  Timoshenko  beam  theory. 
Beam  motions  are  governed  by  the  following  equations  [lief.  4|. 

a2  „  ,  -.2 

a\  (-T-*       =      2  ,'2•2':, 

ex         tx  ct 

,2,  i  «2  i 

2    C^O      .       2,2/    °y  ,,         C   0  ,-,   ,0v 

a2 — r  +  ^k0{-^-(f))  =  — r-  (2.28) 

dx2  GX  dr 

kG*  E*  I 

where  a\  =  — - — ,  a\  —  ——  and  kl  —  —  The  complex  moduli  G"  and  E'  are  defined  as: 

G*  =  G(1  +uh.,)  (2.29) 

E*  =  E{1  +  /;/,,)  (2.30) 

where  subscript  vt  is  an  abbreviation  of  Viscoelastic  beam  using  Timoshenko  beam  the- 
ory, and 

G  :  shear  modulus  of  the  beam 

E  ;  elastic  modulus  of  the  beam 

>/lf  ;  loss  factor  of  the  viscoelastic  material 


A  :  cross  sectional  area  of  the  beam 
/  :  moment  of  the  inertia  of  the  cross  sectional  area 
p  :  density  of  the  beam  material 
k  :  shear  deformation  factor 

The  general  forms  y(x,t)  and  <f)(x,t)  for  the  propagating  harmonic  wave  in  a  Ti- 
moshenko  beam  are 

y(x,r)  =  YQeKk*x-Mt)  (2.31) 

0(jc,O  =  %ei{kXx-a,I)  (2.32) 

where  k"  is  the  complex  wave  number,  and  co  is  the  frequency  of  the  propagating  wave. 
Inserting  the  wave  equations  (2.31  and  2.32)  into  the  governing  beam  equations  (2.27 
and  2.28)  becomes 

(co2  -  a{(k*)2)  Y0  +  (  -ia2kxy%  =  0  (2.33) 

{ia2k%k*)Y0  +  {co2  -  a2k^  -  4(k*)2)%  =  0  (2.34) 

The  dispersion  relation  of  waves  is  obtained  from  equations  (2.33  and  2.34) 

a[a2{kx)   —  (a{  +  a2)co~(k*)   +  (co   —  alk0co  )  =  0  (2.35) 

with  the  relationship  between  Y0  and  cl>0. 

pu)   -kG  (k   ) 

°o  = ~ }  o  =  #0 }  o  (•-•  -6) 

ikG  k 

where 

pgj2  -KG*(kx)2 

^o  = ~ —  ^•J') 

ikG  k 

The  equation  (2.35)  gives  the  four  different  wave  numbers  as  function  of  wave  frequen- 
cies. 


[a2  +  a2)<oJ  +  ^  (rtf  +  a2):o)4  -  afo^w4  -  a^co2) 
fc,  =  J * -j-2 (2.38.a) 

2fl]a2 


2  2\     2  i  2\?     4  2^4  ^i-1 

(a,  +  a*,  )co  +  \/  (a,  +  a-, )  to   —  a,  a'{(0   —  a,  knCJ  ) 

A2X  =  -    /  — = ^ — ■ ^ LJ _! (2.3S.Z?) 

V  2a7a-, 


1/2  2\     2  2  '  2     4  ''2-1  2    2^ 

/  (a,  +  a-))(o  —  x  (a,  +  cfo )  co   —  a,a-.(co   —  a,  /c(1co"i 

A3X  =     /— ^— i =^-r — — (2.38.c) 


•2  2       2  '  /    2  2s  2     4  2    2/4  2    7     2- 

( 67,' +  c7t  )co    —  v'  (rtf  +  a2)  CO    —  c7J'a;(to    —  aJ"/c,^o") 
A4  =  -  a  /  -  — —  "  (2.38.d) 

Therefore,  the  equations  (2.31  and  2.32)  are  expanded  using  kf,  A:\  A;  and  A; 

y(xj)  =  ( }>/A:>  "x  +  Y2eik*x  +  Y,elk*x  +  YA^x)e-'l<al  (2.39) 

4>{x.i)  =  ((D^^  +  ®2eik*x  +  <&2elk*x  +  ®±eik*x)e-i(at  (2.40) 

The  equation  (2.40)  becomes  the  following  equation 

O(x)  =  Rj  F, A*  +  K2 h^X  +  R3  Y3etk*x  +  #4 *4«^<JC  (2.4 1 ) 

where 


,~X,  ,_X     ■ 


i        /■-*■,/> 
pco   —  k(j  (A., , 

/?„  = -^-  (2.42) 

/K.'(j     K/2 

for«=  1.  2,  3.  4. 

For  a  finite  viscoelastic  beam  with  a  length  /,  which  is  free  at  x  =  I  and  excited 
at  x  =  0,  four  boundary  conditions  are  defined  as: 

y(0)=y  (2.43) 

0(0)  =  0  (2.44) 

v    86(1) 
:V/(/)  =  £:x/-tJ~--0  (2.45) 

a 

S([)  =  kG*A( -~-<f>)  =  0  (2.46) 

From  those  boundary  conditions,  the  constants  }',,  Y2,  l\  and  )\  are  determined 
as    functions    of  y    and    </>.    At    the    driving    point,    x  =  0.    input    transverse    force, 
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F(0)  =  -S(0),  can  be  calculated  using  equation  (2.46),  which  can  be  expressed  using 
impedances  and  velocity  at  the  driving  point  as  follows: 

F(0)  =  Z,  V{0)  +  Z2O(0)  (2.47) 

where  the  complex  transverse  velocity.  V,  and  the  complex  angular  velocity.  <L>,  are  re- 
lated to  the  corresponding  displacement  as  follows: 

J-  =  ±  =  -iM  (2.48) 

y       (p 

If  there  is  no  angular  displacement,  <f>,  at  the  center  of  the  beam,  that  is,  0(0)  =  0,  the 
equation  (2.47)  becomes 

F(0)  =  Z]V(0)  (2.49) 

Therefore,  the  impedance  of  the  Timoshenko  beam  at  the  driving  point,  Z/w,  can  be  cal- 
culated using  the  following  equation. 

F/Vr(0)     1 

where  subscript  fit  is  an  abbreviation  of  the  Finite  Viscoelastic  beam  using  Timoshenko 
beam  theory. 

3.     FINITE  CONSTRAINED  LAYER  BEAM 

In  this  section,  impedances  of  the  constrained  layer  beam,  which  has  viscoelastic 
material  between  the  two  elastic  layers,  are  studied.  The  theory  of  a  damped  sandwich 
beam  was  developed  by  R.A.  Ditaranto  [Ref  5],  who  extended  Kenvin's  [Ref.  6  ]  bas- 
ically similar  analysis.  D.J.  Mead  and  S.  Markus  [Ref.  7]  expanded  the  concept  by  con- 
sidering these  earlier  works.  The  result  of  their  elforts  is  a  sixth  order  differential 
equation  of  motion  which  is  expressed  in  terms  of  the  transverse  displacement,  y,  and 
longitudinal  displacement,  u.  for  the  constrained  layer  beam.  The  constrained  layer 
beam  will  be  excited  harmonically  at  its  center  by  a  shear  force.  S,  as  shown  in 
Figure  3  on  page  12.  This  equation  of  motion  leads  to  the  impedances  of  the  sandwich 
beams  with  six  boundry  conditions,  which  are  used  for  evaluating  impedances  depending 
on  frequency,  viscoelastic  material  properity,  and  two  elastic  material  properties. 

From  the  [Ref.  7],  the  governing  sixth  order  differential  equation  is 
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9 


constrained  laver 


w////////  v^scoe^as^c  material y/j//////. 


constrained  layer 


Figure  3.      Coordinates  and  sign  conventions  for  a  semi-finite  constrained  laver  beam. 


~6  -,4  ,2 

V    V  C    V  1  0  p 

T-r-^(i  +  Z)-rT  =  -^-(T-V-z)^) 

CX  CX  ud        OX 


(2.51) 


where 


r   v 
p  =  -m  -fT 
ct 


(2.52) 


with  the  relationship  between  y{xj)  and  u{xj)  expressed  as: 


e2tt       0  BuDj    8y 

-7-7-  -  Bbu  =  - 


ox 


E3H3CC    ex 


(2.53) 


The  parameters  in  this  equation  are  defined  as  follows: 


b    '  H2  {  E[H[       E3H3  ' 


[2.54.0) 
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Z  =  — ( 1  (2.54./>) 

Z?d  l  £,//,  + £3//3  J  7 

Dd=ElIl  +  E3I3  (2.54.0 

Q  =  //2  +  ^  (//,+//,)  (2.54*0 

/,-- j£-  (2.54.,) 

/3  =  -~-  (2.54./) 

G*  =  G(1  +/;?c)  (2.55) 

where  subscript  c  is  an  abbreviation  of  constrained  layer  beam,  and 

0  :  shear  modulus  of  the  viscoelastic  layer 

E  :  elastic  modulus  of  the  elastic  layer 

)]c  :  loss  factor  of  the  viscoelastic  layer 

/  :  moment  of  the  inertia  of  the  cross  sectional  area 

H  :  thickness  of  the  constrained  layer  beam 

m  :  mass  per  unit  length  of  the  whole  three  layer  section. 

The  general  form  of  transverse  displacement,  y(x,t),  and  longitudinal  displace- 
ment, u(x,t),  for  the  propagating  harmonic  waves  in  a  constrained  layer  beam  are 

y{x.t)=  V  (2.56) 

u{x,t)  =  UQe  (2.;w) 

where  k*  is  a  complex  wave  number  and  to  is  the  frequency  of  the  propagating  wave. 
After  inserting  the  wave  equations  (2.56  and  2.57)  into  the  sixth  order  beam  equation 
(2.51),  the  dispersion  relation  is  obtained  as  follows: 

(k*  f  +  Bb(\  +  Z)(k*f  -  -22-  (kxf  -  J^_  =  0  (2.58) 
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There  are  six  different  roots  with  six  different  complex  wave  numbers  in  the  above 
equation,  which  means  a  wave  of  six  different  modes  propagate  through  the  constrainei 
layer  beam.  From  the  equations  (2.56  and  2.57).  we  may  rearrange  equation  (2.53) 

U0  =  R0Y0  U.59) 

where 

iBbZDdk» 
R0  = ^—^ (2.60) 

E3H3Cc((k  ] I  +  Bb) 

The  general  form  of  the  transverse  displacement  equation  (2.56),  y,  and  longi- 
tudinal displacement  equation  (2.57),  u  ,  are  obtained  using  the  six  complex  wave  num- 
bers and  a  ratio  R0.  The  equations  (2.56.  2.57  and  2.59)  as  follow: 

;      ,\        cv     lKx    ,     v     IkiX   ,     v     ik$x   .     v-     <ki-x   ,     V     ^<x    <     v     &6X\    —loot  ,  -,  s  ,  . 

y{x,t)  =  {lle  '    -f-  i2e  -    +  Y3e  3    -I-  1+e  4    +  1 5e  -    +  i6e  *  )e  (_.6i) 

u{xj)  =  (RlYleik'x+R2Y2eik'x+R2Y3ek'x+RAY^x 

+  R5Y5elk*x  +  R6Y6eikhe-ia}l  (2.62) 

where 

iBbZDdk* 

K  =  —-~ — (2-63) 

((A,*)2  +Bb)E3H3Cc 

n—  1,  2.  3.  4.  5.  6  For  a  constrained  layer  beam  with  a  finite  length  /  which  is  free 
at  x  =  /  and  excited  at  one  end  at  x  =  0,  the  six  boundary  conditions  are  defined  as  fol- 
lows [Ref.  7|: 

y{Q)=y  (2.64) 

rfvfO') 


ux 


=  6  (2.65) 


m(0)  =  w  (-.66) 


c3y       „  _  8y        DdBbZ      £,//,  +  /:,//, 


5(/)  =  Dd  -rr  -  DdBbZ  -f-  +  -y-  (  \,       "       )u  =  0  (2.67) 


OX  u-v  W  M'M 
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M(l)  =  Dd-^-~  C.E3H3  4^  =  0  (2.6S) 

cx  °  -v 

X(f)=  E3H3-^-  =  i)  (2.09) 

cx 

The  six  constants  1'  (7=1,  2.  3,  4,  5  and  6)  are  determined  from  these  six  boun- 
dary conditions  and  are  a  function  of y,  6  and  u.  At  the  driving  point,  x  =  0,  the  input 
transverse  force.  F(0)  =  —  S(0).  moment.  A/(0),  and  normal  force.  A'(0),  can  be  calculated 
using  equations  (2.67.  2.6S  and  2.69).  These  input  force  and  moment  can  be  expressed 
using  impedances  and  velocities  and  angular  velocity  at  the  driving  point  as  follows: 

F(0)  =  Z,  V{0)  +  Z2Q{0)  +  Z3  U{0)  ( 2. 70) 

Where  the  complex  transverse  velocity,  V,  and  the  complex  angular  velocity, 
Q.,  and  the  complex  longitudinal  velocity,  U,  are  related  to  the  corresponding  displace- 
ment as  follows: 

—  =  -r  =  — =  -/co  (2.71) 

If  there  is  no  rotational  displacement,  6,  and  no  longitudinal  displacement,  u.  at  the 
center  of  the  beam,  Zx  represents  the  impedances  of  the  transverse  force  at  the  driving 
point  of  the  constrained  layer  beam.  That  is,  Q(0)  =  0  and  [/(0)  =  0.  Therefore,  the 
equation  (2.70)  becomes 

F(Q)  =  ZlV(0)  (2.72) 

Therefore,  the  impedances  of  the  constrained  layer  beam  at  the  driving  point.  Zc,  can 
be  calculated  using  the  following  equation. 

FJ0)    i 
Zm-^f  (2.73) 

Where  subscript  c  is  an  abbreviation  of  the  constrained  layer  beam. 

C.     COMPUTER  IMPLEMENTATION 

The  computer  programs  using  Reduce  and  FORTRAN  languages  are  used  for  cal- 
culation of  the  theoretical  impedances  from  the  beam  theories  in  this  paper.  Reduce  was 
developed  by  Anthony  C.  Hearn  [Ref.  8],  and  is  a  programming  language  which  solves 
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algebraic  operations  non-numerically.  It  can  manipulate  polynomials  in  a  varity  of  forms 

(substitution,  differentiation  and  integration)  and  solves  one  or  more  simultaneous  al- 
gebraic equations.  After  the  impedances  were  evaluated  by  a  Reduce  program  with 
constant  variables,  a  FORTRAN  program  calculated  the  impedance  with  real  dimen- 
sions and  properties  of  the  viscoelastic  and  constrained  layer  beam  waveguide  absorber 
for  the  0-2000Hz  frequency  range.  All  constant  variables  of  the  Reduce  program  were 
defined  in  each  FORTRAN  program.  We  tried  to  use  the  same  notations  for  the  Reduce 
program,  the  FORTRAN  program  and  the  text  (Appendix  A). 

Three  Reduce  programs  and  a  FORTRAN  program  were  used  to  study  each  theory. 
The  first  Reduce  program  defined  the  moment,  force  and  displacements  which  included 
the  constants  )".,  where  n  equals  1  and  3  for  an  infinite  viscoelastic  beam;  n  equals  1.  2, 
3  and  1  for  a  finite  viscoelastic  beam  and  n  equals  1,  2.  3.  4.  5.  and  6  for  a  constrained 
layer  beam.  Yn  can  be  solved  with  boundary  conditions.  After  running  the  first  Reduce 
program,  new  constant  terms  are  defined  for  Y„.  The  second  Reduce  program  includes 
those  new  constant  variables  which  were  the  results  of  the  first  Reduce  program.  The 
second  Reduce  program  gives  the  constant  variables  Yn.  Finally,  the  third  P^educe  pro- 
gram can  evaluate  the  impedance  of  the  beam  with  longitudinal  and  angular  displace- 
ment (and  normal  displacement  for  the  constrained  layer  beam).  The  process  o[  the 
constrained  layer  beam  is  similar  up  to  the  second  Reduce  program.  However,  the 
computer  finds  it  hard  to  handle  6  simultaneous  equations  with  6  non-numerical  vari- 
ables, Y„.  Therefore,  we  solved  }"4,  Y5  and  Y6  in  the  second  Reduce  program  and  then 
substitute  these  values  to  }',,  Y2  and  Y3  in  the  third  Reduce  program. 

Wave  numbers  as  functions  of  frequencies,  i.e.  dispersion  relation,  of  waves  in  wav- 
eguide absorbers  were  described  in  previous  sections.  It  was  easy  to  find  the  wave  num- 
ber of  a  viscoelastic  beam  through  Bernoulli-Euler  beam  theory  and  Timoshcnko  beam 
theory  in  equations  (2.12  and  2.38).  Flowever.  wave  numbers  of  the  constrained  layer 
beam  were  expressed  with  a  complex  six  order  polynomial.  Therefore,  they  were  calcu- 
lated using  Newton's  method.  The  equation  (2.58)  is  rearranged  using  real  constants 
/;,.  cr  and  dr  and  imaginary  constants  /?,  and  J,. 

{K*f  +  {br  +  ibi){Kx)2  +  cr(Kx)  +  (dr  +  id))  =  0  (2.74) 

where 

^x  =  (k*)2  (2.15.0) 
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br  =  Re(Bb  +  BbZ)  (2.75.*) 

/,.=  Im(Bb  +  BbZ)  (2.7 5.  c) 

i 
cV  =  -^—  (2./5.«/) 

dr  =  Re{ L )  (2.7:>.<0 

—Biinto 

di  =  Im( - )  (2.75/) 

Let  the  complex  variable  K"  =  t  +  iq  [t  and  q  are  real  numbers).  The  equation  (2.7-1) 
becomes  two  equations  coupled  with  two  variables  (t  and  q). 

J[t,q)  =  t 3  +  b,.t2  +  (  -3<72  -  2%  +  cr)r  +  (  -brq2  +  <afr)  ( 2.76) 

g(t,q)  =  ./3  +  biq2  +  (  -3r2  -  lbrt  -  c^  +  (  -L^t2  -  d-J  (2.77) 

Starting  from  initial  guesses  (t0,  qn)  which  are  near  the  roots  iteration  was  continued  until 
either  the  successive  t  and  q  values  converged  to  the  certain  values  or  the  value  of  the 
function  is  sufficiently  near  zero.  After  a  few  iterations  of  equations  (2.76  and  2.77) 
using  a  FORTRAN  program  (Appendix  A),  we  can  get  three  different  square  of  wave 
number  values.  K\  were  obtained.  So,  finally,  this  method  calculates  six  different  wave 
numbers,  k;,  from  +  v  K*  ,  where  n  equals  I,  2,  3,  4.  5,  and  6. 

Impedances  of  a  20"  viscoelastic  beam  were  calculated  using  these  computer  pro- 
grams for  viscoelastic  beam  properties  of  200,000  psi  constant  shear  modulus  and  0.2 
energy  loss  factor.  The  results  of  these  simulations  are  shown  in  Figure  -1  on  page  19 
and  Figure  5  on  page  20  which  reflect  sharp  resonance  frequencies.  These  sharp  reso- 
nance phenomena  reduced  for  the  same  beam  with  higher  energy  loss  factor  of  0.5.  For 
cases  using  the  Bernoulli-Euler  beam  theory,  the  impedances  of  the  finite  viscoelastic 
beam  approached  to  the  infinite  viscoelastic  beam  waveguide  absorber  with  increasing 
energy  loss  factor  and  with  increasing  frequency.  The  impedance  using  Timoshenko 
beam  theory  differs  from  the  impedance  from  the  Bernoulli-Euler  beam  theory 
(Figure  6  on  page  21  and  Figure  7  on  page  22),  because  Timoshenko  beam  theory 
considers    shear    deformation    and    rotarv    inertia    effect    of   the    viscoelastic    beam. 
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Impedances  of  a  20"  constrained  layer  beam  with  viscoelastic  layer  properties  of  200  psi 
constant  shear  modulus  and  0.5  and  1.0  energy  loss  factors  were  calculated.  Figure  S 
on  page  23  and  Figure  9  on  page  24  compare  the  impedances  vs.  frequency  depending 
on  the  energy  loss  factor  of  the  constrained  layer  beam.  The  impedance  of  the  high 
damping  constrained  layer  beam  shows  the  smooth  and  shifted  resonance  frequencies 
compared  with  the  low  damping  constrained  layer  beam. 
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Figure  4.      Real  part  of  the  16"  viscoelastie  beam  waveguide  absorber  theoretical  im- 


pedances (r/  =  0.2) 
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Figure  5.      Imaginary  part  of  the  16"  viscoelastic  beam  waveguide  absorber  theore- 


tical impedances  (;/  =  0.2)  . 


20 


u 
z 

D 
CT 

CO 
> 

a 
o 

z 

< 

a 
w 
cu 


3'T  0*1  9'0 

(Nl/OaS-ifll)  'CKMWI  1V3H 


Figure  6.      Real  part  of  the  16"  viscoelastic  beam  waveguide  absorber  theoretical  im- 
pedances (?/  =  0.5)  . 
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Figure  7.      Imaginary  part  of  the  16"  viscoelastic  beam  waveguide  absorber  theore- 
tical impedances  (rj  =  0.5)  . 
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Figure  8.      Real  part  of  the  20"  constrained  layer  beam  waveguide  absorber  theore- 
tical impedances. 
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Figure  9.       Imaginary  part  of  the  2U"  constrained  layer  beam  waveguide  absorber 
theoretical  impedances. 
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III.     EXPERlMEiNT 

A.     IMPEDANCES  OF  THE  TEST  PLATE 

A  rectangular  aluminum  plate  with  clamped  boundary  condition  was  selected  as  the 
test  structure,  because  it  has  enough  closely  space  modes  in  the  frequency  range 
(100Hz-2Kz).  An  aluminum  block  (30"  x  22"  x  2')  was  carved  inside  to  make  a  24"  x  16" 
x  0.3125"  plate  with  clamp  boundary  condition. 

The  impedances  of  the  test  plate  were  measured  from  impact  hammer  tests.  The 
aluminum  plate  was  excited  by  a  PCB  0S6B03  impact  hammer.  Input  impact  force  was 
measured  by  a  force  transducer  at  the  tip  of  the  impact  hammer.  Responses  were  meas- 
ured by  an  Endevco  225OA-10  piezoelectric  accelerometer  which  was  attached  under  the 
impact  point  for  locations  1,  2  and  3  (Figure  10  on  page  26  and  Figure  11  on  page 
27).  Location  1  is  the  center  of  the  plate.  Location  2  is  the  symmetric  point  of  location 
4,  where  a  shaker  will  be  attached  during  the  damping  measurement  tests.  Location  3  is 
a  general  point  which  is  not  on  any  line  of  symmetry.  A  baseband  measurement  for  each 
location  was  taken  from  0  to  2000Hz.  An  average  of  15  data  samples  was  taken  for  each 
measurement  and  analyzed  using  a  HP-3562A  dynamic  signal  analyzer.  Impedance  of 
the  aluminum  plate  at  each  location  was  determined  after  a  predetermined  data 
processing  procedure  using  the  HP-3562A  dynamic  signal  analyzer. 
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Figure  10. 


Arrangement  for  test  plate  impedance  measurement. 


26 


Figure   11.      Test  plate  impedance  measurement  configuration. 
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B.     IMPEDANCES  OF  THE  WAVEGUIDE  ABSORBER 

Two  viscoelastic  beams  and  two  constrained  layer  beams  were  selected  as  the  wav- 
eguide absorber  specimens.  Dimensions  and  physical  properties  are  shown  in  Table  1. 
Mechanical  properties  of  a  viscoelastic  material  depend  on  temperature  and  vibration 
frequencies.  The  loss  factors  and  shear  modulli  of  the  viscoelastic  material  used  in  the 
specimens  at  room  temperature  (15°F)  are  shown  in  Figure  12  on  page  30  and 
Figure  13  on  page  31.  These  data  were  obtained  from  Design  Data  Sheets  in  a  book  by 
Nashif,  et.  al.  [Ref.  9].  Analytic  expressions  of  the  loss  factors  and  shear  moduli  were 
determined  from  curve  fits  to  these  data.  The  following  four  equations  represent  the  loss 
factors  and  shear  moduli  of  the  LD-400  and  ISD-112  viscoelastic  material. 

,/v  =  0.65,[-°-52732|los^)l1'"1  .3.1) 

^=llt;[-0.3906|log(^)|^,  (32j 

Gv  =  0.00002503/'  -0.1752/  +457.5883/+29280  (3.3) 

Gc  =  -0.00001258/  +0.2472/ +  74.988  (3.4) 

f  represents  the  frequency.  ?■/„  and  ;/,  represent  the  loss  factor  of  LD-400  and  ISD-112 
viscoelastic  material.  G„  and  Gc  represent  the  shear  modulus  oi"  LD-400  and  ISD-112 
viscoelastic  material.  Figure  14  on  page  32  shows  waveguide  absorbers  which  were  used 
in  this  study. 

Impedances  of  waveguide  absorbers  at  the  attachment  points  were  determined  from 
steady  state  random  vibration  tests.  The  waveguide  absorber  was  excited  by  a  Wilcoxon 
F4,  F7  vibration  generator  using  a  band  limited  white  noise  signal  (0-2000Hz)  from  a 
HP-3562A  dynamic  signal  analyzer  (Figure  15  on  page  33).  The  waveguide  absorber  was 
mounted  on  the  Wilcoxon  F4,  F7  vibration  generator  (Figure  16  on  page  34).  Force  and 
acceleration  were  measured  using  a  force  transducer  and  a  piezoelectric  accelerometer  in 
the  impedance  head  of  the  Wilcoxon  F4,  F7  shaker.  These  signals  were  analyzed  by  a 
HP-3562A  dynamic  signal  analyzer.  100  data  samples  were  averaged  for  each  measure- 
ment. 

The  measurement  of  impedances  included  the  impedances  of  bolts  and  nuts  and 
connecting  aluminum  pieces.  Therefore,  the  impedances  of  bolts,  nuts,  and  connecting 
aluminum  pieces,  the  so-called  mass  effect,  were  measured  separately.  Total  impedances 
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minus   mass   effect   represented   the   impedances   of  the   waveguide   absorbers.    These 
experimental  results  are  compared  with  the  theoretical  predictions  in  Chapter  IV. 


Table   1.       DIMENSIONS  AND  PROPERTIES  OF  THE  WAVEGUIDE  ABSOR- 
BER SPECIMENS. 


SPECI- 
MEN 1 

SPECI- 
MEN 2 

SPECI- 
MEN 3 

SPECI- 
MEN -4 

WAVEGUIDE  ABSOR- 
BER TYPE 

VI  SCO- 
ELASTIC 

VI  SCO- 
ELASTIC 

CON- 
STRAINED 

LAYER 

CON- 
STRAINED 
LAYER 

MATE- 
RIAL 

VI  SCO- 
ELASTIC 

LD-400 

LD-400 

ISD-112 

ISD-112 

CON- 
STRAINED 

LAYER 

Al 

Al 

LENGTH  (in) 

16.0 

20.0 

16.0 

20.0 

WIDTH  [in) 

1.0 

1.0 

1 .0 

1.0 

THICK- 
NESS 

(in) 

VI  SCO- 
ELASTIC 

0.35 

0.35 

0.065 

0.065 

CON- 
STRAINED 

LAYER 

0.0625 

0.0625 

DEN- 
SITY 

(  lbm  ) 

VI  SCO- 
ELASTIC 

0.0550.5 

0.05505 

0.03663 

0.03663 

CON- 
STRAINED 

LAYER 

0.09832 

0.09832 
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Figure   12.      Loss  factor  vs.  frequency  for  a  LD--4U0  and  a  ISD-112. 
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13.      Shear  modulus  vs.  frequency  for  a  LD-400  and  a  ISD-112. 


31 


-        ;•     .::■-.-..■;.     ■*>,      - 


.  ■'  '■    ■ 


•  4 


Or  » 


^ 


iiis 


Figure   14.      Waveguide  absorbers. 
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Figure   16.      VVilcoxon  F4/F7  shaker  ^vith  waveguide  absorber  mounted. 
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C.     DAMPING  OF  THE  PLATE 

Modal  damping  values  of  the  plate  were  determined  from  structural  transfer  func- 
tions at  a  driving  point  and  between  input  and  output  points  on  the  test  plate  which 
were  measured  using  random  excitations.  The  HP-3562A  dynamic  signal  analyzer  was 
used  in  the  measurement  of  the  transfer  functions.  The  Wilcoxon  F3  vibration  generator 
excited  the  plate  at  location  4  using  random  oscillation  signals  from  a  HP-3562A  dy- 
namic signal  analyzer.  Input  forces  and  output  responses  at  the  driving  point  were 
measured  using  a  force  transducer  and  a  piezoelectric  accelerometer  in  the  impedances 
head  of  Wilcoxon  F3  shaker  (Figure  17  on  page  36).  The  damping  measurement  con- 
figuration is  shown  in  Figure  IS  on  page  37.  It  was  necessary  to  measure  responses  at 
several  points  other  than  the  driving  point  to  accurately  measure  damping  values  of 
closely  spaced  modes.  These  were  measured  using  the  Endevco  2250A-10  piezoelectric 
accelerometer.  Modal  damping  values  of  all  the  modes  in  the  frequency  range  of 
100-2000Hz  were  determined  from  the  response  functions  using  the  Half  power  Band 
Width  method  [Ref.  10]  and  a  curvefit  method  [Ref.  1 1] 

Modal  damping  values  of  the  plate  without  a  waveguide  absorber  were  measured 
as  the  baseline.  A  waveguide  absorber  was  attached  to  the  plate  at  location  1.  The 
waveguide  absorber  was  changed  after  each  measurement  of  the  damping  for  four  dif- 
ferent waveguide  absorbers  which  were  described  in  section  4-B.  After  the  damping  was 
measured  at  location  1,  measurements  were  taken  at  location  2  and  location  3. 

The  most  effective  waveguide  absorber  at  each  location  was  selected.  The  20"  vis- 
coelastic  beam  waveguide  absorber,  the  16"  constrained  layer  beam  waveguide  absorber 
and  the  16"  viscoelastic  beam  waveguide  absorber  were  most  effective  at  location  1,  2 
and  3.  respectively.  Therefore,  the  damping  with  two  waveguide  absorbers  attached  at 
locations  1  and  2,  and  the  damping  with  three  waveguide  absorbers  at  locations  I.  2  and 
3  were  also  measure,  to  see  the  effect  of  multiple  waveguide  absorbers. 
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Figure  17.      Arrangement  for  test  plate  damping  measurement. 
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Figure  18.      Damping  measurement  configuration. 
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IV.     RESULTS  AND  DISCUSSIONS 

A.     IMPEDANCES  OF  THE  TEST  PLATE 

Impedances  of  the  test  plate  were  measured  from  input  hammer  tests  as  described 
in  Section  3. A.  Since  the  aluminum  test  plate  is  a  lightly  damped  structure,  the  impe- 
dances at  resonance  frequencies  are  very  low  and  at  anti-resonance  frequencies  are  very 
high.  Phases  with  ISO0  phase  shift  are  near  ±90°  both  at  the  resonance  and  at  the 
anti-resonance  frequencies. 

Since  location  1  is  the  center  point  of  the  plate,  only  small  numbers  of  resonance 
and  anti-resonance  frequencies  appear  (Figure  19  on  page  39  and  Figure  20  on  page 
40).  Neither  location  2  nor  location  3  are  on  any  line  of  symmetry  of  the  test  plate:  all 
the  resonance  and  anti-resonance  points  appear  in  the  impedance  vs.  frequency  dia- 
grams. (Figure  21  on  page  41.  Figure  22  on  page  42.  Figure  23  on  page  43  and 
Figure  24  on  page  44).  However,  the  anti-resonance  frequencies  at  location  2  and  those 
at  location  3  are  different. 
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Figure  19. 


Impedance  magnitude  of  test  plate  at  location  I, 
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Figure  2U.      Impedance  phase  of  test  plate  at  location  l. 
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Figure  21. 


Impedance  magnitude  of  test  plate  at  location  2. 
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Figure  22. 


Impedance  phase  of  test  plate  at  location  2. 
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Figure  23. 


Impedance  magnitude  of  test  plate  at  location  3. 
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Figure  24.      Impedance  phase  of  test  plate  at  location  3. 
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B.     IMPEDANCES  OF  THE  WAVEGUIDE  ABSORBER 

Two  different  kinds  of  waveguide  absorbers,  a  viscoelastic  beam  and  a  constrained 
laver  beam,  are  investigated  in  this  study.  In  a  previous  study  by  Ungar  and  Williams 
[Ref.  3].  impedances  were  predicted  only  for  the  infinite  viscoelastic  beam  using  Ber- 
noulli-Euler beam  theory  and  the  theoretical  prediction  showed  wide  difference  from  the 
experimental  results.  In  this  study,  this  previous  work  is  extended  to  include  the  impe- 
dances for  the  finite  viscoelastic  beam.  Since  Bemoulli-Fuler  beam  theory  does  not  in- 
clude the  shear  deformation  and  rotary  inertia  effect  of  the  viscoelastic  beam,  the 
predictions  from  the  Timoshenko  beam  theory  are  studied  to  investigate  these  effects  on 
wave  propagation  and  impedances,  and  are  compared  with  the  results  from  the  Ber- 
noulli-Euler beam  theory.  The  impedances  of  the  constrained  layer  beam  waveguide  ab- 
sorber were  predicted  using  the  sixth  order  beam  theory. 

The  impedances  from  these  theoretical  predictions  are  compared  with  the  impe- 
dances from  the  experimental  measurements  for  the  16"  and  20"  length  waveguide  ab- 
sorbers of  two  different  kinds.  In  the  calculation  of  theoretical  predictions,  the  loss 
factors  and  shear  moduli  of  viscoelastic  material  were  varied  as  function  of  frequency 
as  described  in  Section  3.B.  During  the  experiment,  the  Wilcoxon  F4  F7  vibration  gen- 
erator produced  characteristic  errors  which  occured  at  different  frequencies  depending 
on  the  mass  of  the  waveguide  absorbers.  Therefore  the  viscoelastic  beam  has  unreliable 
data  around  700Hz  and  the  constrained  layer  beam  has  unreliable  data  around  3.501  Iz. 

Figure  25  on  page  47  and  Figure  26  on  page  -IS  show  the  three  theoretical  impe- 
dances compared  with  experimental  data.  The  Bernoulli-Euler  beam  theory  for  finite 
beams  shows  some  standing  wave  effect  at  low  frequency  range  due  to  reflection  at  the 
free  end  o[  the  beam,  but  approaches  the  results  of  the  infinite  beam  as  frequency  in- 
creases. The  results  of  the  Bernoulli-Euler  beam  theory  do  not  quite  correspond  to  ex- 
perimental results  for  the  viscoelastic  beam.  Impedances  from  the  Timoshenko  beam 
theory  closely  follow  the  experimental  results  for  viscoelastic  beams  as  shown  in 
Figure  25  on  page  47  and  Figure  26  on  page  48.  These  show  that  the  shear  deformation 
and  rotary  inertia  effects  are  significant  for  the  impedance  prediction  of  beam  type 
waveguide  absorbers. 

Figure  27  on  page  49  and  Figure  28  on  page  50  represent  the  real  and  imaginary 
parts  of  the  16"  constrained  layer  beam  waveguide  absorber  and  Figure  29  on  page  51 
and  Figure  30  on  page  52  represent  the  real  and  imaginary  parts  of  the  20"  constrained 
layer  beam  waveguide  absorber.  In  both  beams,  the  theoretical  results  of  the  impedances 
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using  the  sixth  order  beam  theory  and  the  resonance  frequency  patterns  are  almost  the 
same  as  the  experimental  result. 

Both  Timoshenko  beam  theory  and  the  sixth  order  beam  theory  predict  resonance 
frequencies  shifted  amount  from  the  experimental  results  as  frequency  increases,  and  the 
impedance  magnitudes  of  experimental  results  are  bigger  than  the  theoretical  results. 
These  difference  between  theoretical  and  experimental  impedances  may  derive  from  the 
inaccuracy  of  shear  modulus  and  loss  factors  which  values  were  selected  from  15CF  room 
temperature  and  changed  depending  on  frequency.  However,  the  theoretical  prediction 
using  Timoshenko  beam  theory  for  the  viscoelastic  beam  and  the  sixth  order  beam  the- 
ory for  the  constrained  layer  beam  follow  the  same  trend  (Figure  12  on  page  30  and 
Figure  13  on  page  31 ). 

Impedances  of  waveguide  absorbers  are  very  low  and  do  not  show  any  sharp  vari- 
ation for  all  frequencies  in  the  range  tested,  which  are  quite  different  from  the  impedance 
of  the  test  plate  of  very  lightly  damped  structures. 
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Figure  25.       Real  part  of  the  16"  viscoelastic  beam  waveguide  absorber  impedances 
at  the  center  of  the  beam. 
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Figure  26. 


Imaginary  part  of  the  16"  viscoelastic  beam  waveguide  absorber  impe- 
dances at  the  center  of  beam. 
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Figure  27.       Real  part  of  the  16"  constrained  layer  beam  waveguide  absorber  impe- 
dances at  the  center  of  the  beam. 
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Figure  28.       Imaginary  part  of  the  16"  constrained  layer  beam  waveguide  absorber 
impedances  at  the  center  of  beam. 


50 


Figure  29.      Real  part  of  the  20"  constrained  layer  beam  waveguide  absorber  impe- 
dances at  the  center  of  the  beam. 
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Figure  30.       Imaginary  part  of  the  20"  constrained  layer  beam  waveguide  absorber 
impedances  at  the  center  of  beam. 
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C.     DAMPING  OF  THE  PLATE 

The  purpose  of  this  experiment  is  to  investigate  the  effect  of  light  weight  waveguide 
absorbers  on  damping  increase  in  plate  type  structures  over  a  wide  frequency  range. 
Damping  of  the  test  plate  structure  was  measured  with  an  attachment  of  each  one  of 
four  different  waveguide  absorbers  (16"  or  20"  viscoelastic  or  constrained  layer  beam) 
at  three  different  locations.  The  frequency  response  and  the  modal  damping  values  of  the 
aluminum  test  plate  with  a  waveguide  absorber  were  compared  with  those  of  the  same 
plate  without  any  waveguide  absorber  over  a  wide  frequency  range  (0-2000Hz). 

The  modal  damping  values  of  the  test  plate  without  a  waveguide  absorber  are  shown 
in  Figure  31  on  page  56  and  Figure  32  on  page  57.  They  show  that  more  than  20  lightly 
damped  modes  are  in  the  frequency  range  between  100Hz  to  2000Hz  with  the  lowest 
frequency  195  Hz. 

The  damping  contribution  of  the  same  waveguide  absorber  to  a  same  plate  was  dif- 
ferent depending  on  the  location  of  the  waveguide  absorber  attachment  point.  For  the 
20"  viscoelastic  beam  waveguide  absorber,  there  was  marked  damping  increase  for 
modes  lower  than  1 100Hz  when  it  was  attached  at  location  1  (Figure  33  on  page  58  and 
Figure  34  on  page  59).  However,  it  was  more  effective  for  modes  higher  than  1200Hz 
when  it  was  attached  at  location  2  (Figure  35  on  page  60  and  Figure  36  on  page  61). 
When  it  was  attached  at  location  3  its  contribution  to  plate  damping  increase  was  ef- 
fective for  wider  frequency  range  of  250Hz-1500FIz  (Figure  37  on  page  62  and 
Figure  38  on  page  63).  This  trend  applies  to  other  waveguide  absorbers,  though  the 
magnitude  of  damping  contribution  and  the  frequency  ranges  were  different  for  each 
case  (Figures  in  Appendix  B). 

Different  modal  damping  effects  with  a  different  waveguide  absorber  at  the  same 
location  were  also  observed.  At  location  1,  the  20"  viscoelastic  beam  waveguide  absorber 
was  most  effective  for  the  frequency  range  lower  than  1100Hz  (Figure  33  on  page  5S 
and  Figure  34  on  page  59)  and  the  16"  constrained  layer  beam  waveguide  absorber  was 
most  effective  for  the  frequency  range  higher  than  1200Hz  (Figure  39  on  page  64  and 
Figure  40  on  page  65).  At  location  2,  the  20"  and  16"  viscoelastic  beam  waveguide  ab- 
sorbers showed  moderate  contribution  to  plate  damping  and  the  20"  constrained  layer 
beam  waveguide  absorber  contributed  the  least  amount  to  plate  damping.  For  frequency 
ranges  over  1200Hz,  the  16"  constrained  layer  beam  waveguide  absorber  produced  the 
largest  damping  increase.  At  location  3,  the  20"  and  16"  viscoelastic  beam  waveguide 
absorbers  and  the  20"  constrained  layer  beam  waveguide  absorber  produced  moderate 
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damping  increase  over  a  wide  frequency  range  of  250Hz-1500Hz.  However,  the  16" 
constrained  layer  beam  waveguide  absorber  showed  the  least  damping  increase. 

After  the  selection  of  optimum  damping  at  each  location,  two  wave  guide  absorbers 
(20"  viscoelastic  and  16"  constrained  layer  beam)  were  attached  at  locations  I  and  2  and 
three  waveguide  absorbers  (20"  viscoelastic,  16"  constrained  layer  and  16'  viscoelastic 
beam)  were  attached  at  locations  1,  2  and  3  to  see  the  effect  of  multiple  waveguide  ab- 
sorbers on  the  test  plate  damping  (Figure  41  on  page  66).  Figure  42  on  page  67  and 
Figure  43  on  page  68  represent  frequency  response  and  modal  damping  of  the  plate  with 
two  waveguide  absorbers.  Figure  44  on  page  69  and  Figure  45  on  page  70  represent 
frequency  response  and  modal  damping  of  the  plate  with  three  waveguide  absorbers. 
Damping  of  the  plate  increased  with  increasing  number  of  the  waveguide  absorbers 
(Figure  42  on  page  67  and  Figure  44  on  page  69).  The  damping  increase  resulting  from 
attachment  of  two  or  more  absorbers  resemble  the  sum  of  the  seperate  result  of  each. 

Impedance  of  the  test  plate  at  location  1  is  compared  with  the  driving  point  impe- 
dance of  the  20"  viscoelastic  beam  in  Figure  46  on  page  71.  This  figure  shows  that  at 
resonance  frequencies  the  impedance  of  the  test  plate  is  much  smaller  than  that  of  the 
waveguide  absorber.  Figure  47  on  page  72  and  Figure  4S  on  page  73  show  the  com- 
parsion  between  the  impedances  of  the  test  plate  at  location  2  and  location  3  with  the 
drving  point  of  the  16"  constrained  layer  beam  and  16"  viscoelastic  beam,  respectively. 
These  figures  show  that  at  low  frequency  range  (<400Hz)  the  impedances  of  the  test 
plate  were  higher  than  those  of  waveguide  absorbers  and  at  high  frequency  range 
(>400Hz)  impedances  of  waveguide  absorbers  were  higher  at  some  resonance  frequen- 
cies. 

From  equation  (2.7),  the  loss  factor  contribution.  ;/,  of  a  waveguide  absorber  is 
proportional  to  the  following  3  terms: 

Rb  (4.1) 

1 


Zb     2 


(4.2) 


Vs      1 

-p1-!2  (-4.3) 

y  m 
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The  first  term  is  determined  from  waveguide  characteristics  only  and  the  second  term  is 
determined  from  the  impedance  ratio  between  the  test  structure  and  a  waveguide  ab- 
sorber at  the  attachment  point.  However,  the  third  term  depends  on  the  interaction  be- 
tween the  test  structure  and  a  waveguide  absorber  and  the  whole  response  of  the  test 
plate.  If  Z,,  is  much  smaller  than  Z.,  then  Vm  will  remain  constant  before  and  after  the 
attachment  of  a  waveguide  absorber.  Previous  studies  [Ref.  2  and  Ref.  3J  were  based 
on  this  assumption  and  gave  misleading  indications  that  experimental  loss  factors  would 
be  much  higher  than  theoretical  predictions.  In  the  present  study,  as  shown  in 
Figure  46  on  page  71,  Figure  47  on  page  72  and  Figure  48  on  page  73.  the  impedances 
of  the  test  plate  were  similar  in  magnitude  to  impedances  of  waveguide  absorbers  at  re- 
sonance frequencies  and  the  loss  factor  contribution  due  to  waveguide  absorbers  could 
not  be  predicted  using  equation  (2.7). 
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Figure  31.      The  driving  point  frequency  response  of  the  test  plate  without  a  wave- 


guide absorber. 
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Figure  32.       Modal  damping  factors  vs.  frequency  of  the  test  plate  without  a  wave- 
guide absorber. 
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Figure  33.  The  driving  point  frequency  response  of  the  test  plate  with  a  ZQ"  viscoe- 
lastic  beam  waveguide  absorber  (dashed)  at  location  I  and  wichout 
(solid). 
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Figure  34.      Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  20"    vis- 
coelastic  beam  waveguide  absorber  at  location  1  and  without. 
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Figure  35.  The  driving  point  frequency  response  of  the  test  plate  with  a  20"  viseoe- 
lastic  beam  waveguide  absorber  (dashed)  at  location  2  and  without 
(solid). 
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Figure  36.      Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  20"  viscoc 
lastic  beam  waveguide  absorber  at  location  2  and  without. 
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Figure  37.  The  driving  point  frequency  response  of  test  test  plate  with  a  2  0"  vis- 
coelastic  beam  waveguide  absorber  (dashed)  at  location  3  and  without 
(solid). 
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Figure  38.      Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  20"  viscoe- 
lastic  beam  waveguide  absorber  at  location  3  and  without. 
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Figure  39.  The  driving  point  frequency  response  of  the  test  plate  with  a  16"  con- 
strained layer  beam  waveguide  absorber  (dashed)  at  location  1  and 
without  (solid). 
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Figure  -40.      Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  16"    con- 
strained layer  beam  waveguide  absorber  at  location  I  and  without. 
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Figure  41.      The  test  plate  with  three  waveguide  absorbers. 
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Figure  -42.  The  driving  point  frequency  response  of  the  test  plate  with  a  20"  viscoe- 
Iastic  (at  location  1)  and  16"  constrained  layer  (at  location  2)  beam 
waveguide  absorber  (dashed)  and  without  (solid). 
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Figure  -43.  Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  20"  viscoe- 
lastic  (location  I)  and  a  16"  constrained  layer  (location  2)  beam  wave- 
guide absorber  and  without. 
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Figure  44.  The  driving  point  frequency  response  of  the  test  plate  uith  a  20"  viscoe- 
Iastic  (at  location  1),  a  16"  constrained  layer  (at  location  2)  and  a  16" 
viscoelastic  (at  location  3)  beam  waveguide  absorber  (dashed)  and 
without  (solid). 
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Figure  45.  Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  20"  viscoe- 
lastic  (location  1),  a  16"  constrained  layer  (location  2)  and  a  16"  viscoe- 
latic  beam  waveguide  absorber  and  without. 
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Figure  46.      Magnitude  impedances  of  the  test  piate  at  location  1  (solid)  and  of  the 
20"  viscoelastic  beam  at  the  driving  point  (dashed). 
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Figure  47.      Magnitude  impedances  of  the  test  plate  at  location  2  (solid)  and  of  the 
16 ''constrained  layer  beam  at  the  driving  point  (dashed). 
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Figure  48.      Magnitude  impedances  of  the  test  plate  at  location  3  (solid)  and  of  the 
16"  viscoelastic  beam  at  the  driving  point  (dashed). 
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V.     CONCLUSIONS 

The  purpose  of  this  study  was  to  develop  two  kinds  of  high  damping  beam  wave- 
guide absorbers  for  application  to  naval  structures  over  a  wide  frequency  range 
(100Hz- 2000Hz)  and  to  evaluate  theoretical  prediction  schemes  using  various  beam 
theories.  Comparison  between  theoretical  predictions  and  experimental  results  show  that 
Timoshenko  beam  theory  and  the  six  order  beam  theory  can  be  used  in  the  prediction 
of  impedances  of  viscoelastic  beam  waveguide  absorbers  and  constrained  layer  beam 
waveguide  absorbers,  respectively. 

Application  of  waveguide  absorbers  to  a  test  plate  structure  showed  that  damping 
of  the  test  structure  can  be  increased  significantly  with  one  waveguide  absorber  over  a 
frequency  range.  The  magnitude  and  the  frequency  range  of  damping  increase  depend 
on  the  impedance  of  the  waveguide  absorber  and  on  the  location.  The  prediction  of 
damping  increase  due  to  an  attached  waveguide  absorber,  equation  (2.7),  does  not  show 
the  effect  of  waveguide  absorber  completely  since  this  equation  does  not  include  the 
structure  vibration  reduction  effects  due  to  the  waveguide  absorber.  Since  the  contrib- 
ution of  each  waveguide  absorber  to  the  damping  of  the  test  plate  is  additive,  the  modal 
damping  values  of  all  the  vibration  modes  of  the  test  plate  in  a  wide  frequency  range  can 
be  increased  using  multiple  waveguide  absorbers. 
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VI.     RECOMMENDATIONS 

From  the  results  of  present  study  following  are  recommended  for  future  studies  on 
waveguide  absorbers: 

1.  Investigation  of  the  effect  of  a  waveguide  absorber  on  the  reduction  in  vibration 
of  the  original  structure  and  derivation  of  the  relation  between  the  impedances  and  the 
mode  shape  of  the  structure, 

2.  Theoretical  studies  on  impedance  prediction  of  waveguides  of  different  shape  such 
as  circular  viscoelastic  plates  and  application  to  damping  increase, 

3.  Investigation  of  interaction  between  waveguide  absorber  for  multiple  waveguide 
absorbers  applications,  and 

4.  Optimization  of  application  of  waveguide  absorbers  in  design,  number,  and 
location. 
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APPENDIX  A. 

1.     THE  IMPEDANCE  OF  THE  INFINITE  VISCOELASTIC  BEAM  USING 
BERNOULLI-EULER  BEAM  THEORY 


THE  THREE   REDUCE   PROGRAMS   WERE   USED   FOR  EVALUATING  THE    IMPEDANCE   * 

*  OF  THE    INFINITE   VISCOELASTIC   BEAM  USING  BERNOULLI-EULER  BEAM  * 

*  THEORY. 

»'-  y*  .j-  <j~  j-  *»-  V?*V  Vf  ■*•  'V  J-  *V  *V  V?  Vr  V*  Vr  -V  V  r  ~V  ^r  *  V  V?  *V  V?  *V  **-  Vr  ~V  -V  *f*  Vr  ■/■  *'-  ■*•  ■*"  -V  Vr  *V  V*  Vr  Vr  Vr  V?  V~  ■*•  Vr  "V  Vc  *V  ~V  V?  ~V  Vr  **■  iV  Vr  V*  V*  V"  *V  ~ V  "V  Vr  Vf  i'r  #V  V*  -V  V? 

c 

C  THE   FIRST  REDUCE   PROGRAM  DEFINED  TRANSVERSE   FORCE   AND  TOO  BOUNDARY 

C  CONDITIONS. 

C 

YY:  =Y1*EXP( I*K1*X)+Y3*EXP( I*K3*X); 

PP:=DF(YY,X); 

MM:=-EI*DF(PP,X); 

FF:=-DF(MM,X); 

YO:=SUB(X=0,YY)-Y; 

PO:=SUB(X=0,PP)-P; 

BYE; 

C 

C  THE   SECOND  REDUCE   PROGRAM  SOLVED  CONSTANT  VARIABLES    (Yl   AND  Y3). 

C 

Y0:=Y1+Y3-Y; 

P0:=P1*Y1+P3*Y3-P; 

SOLVE(LST(Y0,P0) ,Y1,Y3); 

BYE; 

C 

C  THE  THIRD  REDUCE   PROGRAM  EVALUATED  THE  TRANSVERSE   FORCE. 

C 

YY:  =Y1*EXP( I*K1*X)+Y3*EXP( I*K3*X); 

FFF:=EI*DF(YY,X,3); 

FF:=SUB(X=0,FFF); 

BYE; 


J>  J-  y-  -r^  j*  j„  jm  .r.  j^,  ^f.  tj^  j,  j,  j„  j„  j.  j,  j„  j,  j_  _•„  »r.  j.  .j  „  ^l  J,  *f-fcr*  -•-  J.  „f.  V"'"  ~  '  * '  *'"  *r-  *'*  -,"  "'"  ***  m3~  "**  "'•"  "■'"  *''  V*  *'*■  *'"  ",*  *'-  -'*  "V  **"  V*  "**  *J"  V"  "■'"*  4s  *  V  *'*  "  "V  *V  ■*■ J"  *'*  ~V  *'*  *'* 

*  THIS   FORTRAN   PROGRAM  DEVELOPED  THE    IMPEDANCE   OF  THE    INFINITE 

*  VISCOELASTIC   BEAM  USING  BERNOULLI-EULER  BEAM  THEORY.  * 

DEFINED  VARIABLES; 

COMPLEX  YOUNG'S   MODULUS   OF  THE  VISCOELASTIC   BEAM 

FREQUENCY    (HZ)  * 

TRANSVERSE  FORCE   (LBF) 

COMPLEX   SHEAR  MODULUS   OF  THE   VISCOELASTIC   BEAM 

GRAVITY   (32$    IN/ SEC 

IMPEDANCE   OF  THE   VISCOELASTIC   BEAM  * 

POISSON'S   RATIO 

DENSITY  OF  THE   VISCOELASTIC   BEAM  * 
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J. 

E 

J. 

F 

* 

FF 

* 

G 

tV 

GRAV 

->. 

IMP 

y. 

POISS 

* 

ROH 

*  XB,    XH,XL   :    DIMENSIONS   OF  THE   BEAM   (WIDTH,    HEIGHT,    LENGTH) 
XI  :    AREA  MOMENT  OF   INERTIA 

XK1,    XK2,    XK3,    XK4   :    WAVE   NUMBERS 

*  XNETA        :    ENERGY   LOSS   FACTOR  OF  THE   VISCOELASTIC   BEAM 
VI  :     IMAGINARY   (0,1) 


J*  tlft  J,  »*-  «.'*  iV  iV  »  V  »  V  *tf-  -*f  »' 


frWWWWWWWWWWWftWWWWTPWWWWTVWW^ 


COMPLEX  V1,G,E,KB1,KB2,XK1,XK2,XK3,XK4,P1,P3,A1,A2,A3,A4,FF,IMP 
REAL  F,POISS,GRAV,PI,XB,XH,XG,ROH,XI 
REAL  W , XLOG , XNETA , XA , XCA , XSA , KC , IMPR , IMP I 
C 

OPEN(UNIT=15 ,FILE=' IFI ' ) 
C 

C  CONSTANTS 

C 

F=10.0 
POISS=0.45 
GRAV=386. 0 
PI=3.  1415926 
V1=CMPLX(0.  0,1.0) 
C 

C  DIMENSIONS 

C 

XB=1.  0 
XH=0.  35 
C 

C  PROPERITIES 

C 

ROH=0.05505/GRAV 
XI=(XB*XH**3)/12.  0 
C 

C  FOUR  DIFFERENT  WAVE   NUMBERS 

C 

DO   100    1=1,100 
W=2*PI*F 

XL0G=ABS(L0G10(F/62)) 
XNETA=0.  65*EXP(-0.  52732*XLOG**l.  956) 
C  XNETA=0. 2 

C  XNETA=0. 5 

XG=0.  00002503*F**3-0.  1752*F**2+45  7.  5883*F+29280 
C  XG=200000 

G=XG*CMPLX(  1.  0, XNETA) 

E=2"G*Q+POISS) 

XA=ATAN(  XNETA)/ 4.0 

XCA=COS(-XA) 

XSA=SIN(-XA) 

KC=( SQRT( 1+XNETA**2) )**( -0.  25) 

KB  1=-CSQRT(  ROH*XB*XH*W**2/ ( E*XI ) ) 

KB2=  CSQRT(  ROH*OCB*XH*W**2/(  E*XI ) ) 

XK1=  CSQRT(  KB  1 ) *KC*CMPLX(  XCA ,  XSA) 

XK2=-CSQRT( KB 1)*KC*CMPLX( XCA , XSA) 

XK3=-CSQRT( KB2 )*KC*CMPLX( XCA , XSA) 

XK4=  CSQRT( KB 2 )*KC*CMPLX( XCA , XSA) 
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C     DEFINED  CONSTANTS  FROM  THE  THREE  REDUCE  PROGRAMS 
C 

P1=V1*XK1 

P2=V1*XK2 

P3=V1*XK3 

P4=V1*XK4 

Al=  P3/(P3-P1) 

A2=0 

A3=-P1/(P3-P1) 

A4=0 

FF=-V1*E*XI*( A1*XK1**3+A3*XK3**3 ) 
C 

C     IMPEDANCES 
C 

IMP=2*FF/(V1*W) 

IMPR=REAL( IMP) 

IMPI=AIMAG( IMP) 

WRITE(15,1000)  F,IMPR,IMPI 

F=F+20.  0 
100   CONTINUE 

1000   F0RMAT(1X,D9. 3 ,2X,D12.  5 , 2X,D12.  5) 
CL0SE(UNIT=15) 
STOP 
END 

2.     THE  IMPEDANCE  OF  THE  FINITE  V1SCOELASTIC  BEAM  USING 

BERNOULLI-EULER  BEAM  THEORY 

YY:  =Y1*EXP(  I*K1*X)+Y2*EXP(  I*K2*X)+Y3*EXP(  I*K3*X)+Y4*EXP(  I*K4*X) ; 

PP:=DF(YY,X); 

MM:=-EI---DF(PP,X); 

FF:=-DF(MM,X); 

Y0:=SUB(X=0,YY)-Y; 

P0:=SUB(X=0,PP)-P; 

ML:=SUB(X=L,MM); 

FL:=SU'3(X=L,FF); 

BYE; 

C 

C    THE  SECOND  REDUCE  PROGRAM  SOLVED  CONSTANT  VARIABLES  (Yl,  Y2,  Y3 

C    AND  Y4). 

C 

Y0:=Y1+Y2+Y3+Y4-Y; 

P0:=P1*Y1+P2*Y2+P3*Y3+P4*Y4-P; 

ML:  =M1*Y1+M2*Y2+M3*Y3+M4*Y4; 

FL:  =Fi*Yl+F2*Y2+F3*Y3+F4*Y4; 

SOLVE(LST(Y0,P0,ML,FL),Yl,Y2,Y3,Y4); 

BYE; 

C 

C    THE  THIRD  REDUCE  PROGRAM  EVALUATED  THE  TRANSVERSE  FORCE. 

C 

YY: =Y1*EXP( I*K1*X)+Y2*EXF( I*K2*X)+Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X) ; 

FFF:=EI*DF(YY,X,3); 

FF:=SUB(X=0,FFF); 

BYE; 
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A  * 

*  THIS  FORTRAN  FROGRAM  DEVELOPED  THE  IMPEDANCE  OF  THE  FINITE 
VISCOELASTIC  BEAM  USING  BERNOULLI -EULER  BEAM  THEORY. 

iV  >v 

DEFINED  VARIABLES; 

E       :  COMPLEX  YOUNG'S  MODULUS  OF  THE  VISCOELASTIC  BEAM 

F       :  FREQUENCY  (HZ) 

FF      :  TRANSVERSE  FORCE  (LBF^) 

*  G       :  COMPLEX  SHEAR  MODULUS  OF  THE  VISCOELASTIC  BEAM 

*  GRAV    :  GRAVITY  (323)3  IN/SEC 

IMP     :  IMPEDANCE  OF  THE  VISCOELASTIC  BEAM 

POISS    :  POISSON'S  RATIO 

ROH     :  DENSITY  OF  THE  VISCOELASTIC  BEAM 

*  XB,  XH,XL  :  DIMENSIONS  OF  THE  BEAM  f WIDTH,  HEIGHT,  LENGTH) 
XI      :  AREA  MOMENT  OF  INERTIA 

XK1,  XK2,  XK3,  XK4  :  WAVE  NUMBERS 

XNETA   :  ENERGY  LOSS  FACTOR  OF  THE  VISCOELASTIC  BEAM 

*  VI      :  IMAGINARY  (0,1) 

*  * 

c 

COMPLEX  V1,G,E,KB1,KB2,XK1,XK2,XK3,XK4,P1,P2,P3,P4 

COMPLEX  M1,M2,M3,M4,F1,F2,F3,F4,DEN,A1,A2,A3,A4,FF,IMP 

REAL  F, POISS, GRAV, PI, XB, XL, XH,XG, ROH, XI, W,XL0G, XNETA, XA,XCA,XSA 

REAL  KC,KR1,KR2,KR3,KR4,KI1,KI2,KI3,KI4,IMPR,IMPI 
C 

OPEN( UNIT=15 ,FILE=' FFI ' ) 
C 

C     CONSTANTS 
C 

F=10.0 

POISS=0.45 

GRAV=386.  0 

PI=3.  1415926 

V1=CMPLX(  0.0,1.  0) 
C 

C     DIMENSIONS 
C 

XB=1.  0 

XL=8.  0 

XH=0. 35 
C 

C     PROPERITIES 
C 

ROH=0.  05505/GRAV 
C 

C     FOUR  DIFFERENT  WAVE  NUMBERS 
C 

XI=(XB*XH**3)/12.0 

DO  100  1=1,100 
W=2*PI*F 

XL0G=ABS(L0G10(F/62) ) 
XNETA=0.  65->EXP(-0.  52732 -XL0G **  1.  956) 
C        XNETA=0. 2 
C        XNETA=0.  5 
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XG=0.  00002503*F**3-0. 1752*F**2+457. 5883*F+29280 

C        XG=200000 

G=XG*CMPLX( 1.  0 ,XNETA) 

E=2*G*(1+P0ISS) 

XA=ATAN(XNETA)/4. 0 

XCA=COS(-XA) 

XSA=SIN(-XA) 

KC=(  SQRT(  14XNETA**2))**(  -0.  25 ') 

KB 1=-CSQRT( ROH*XB*XH*W**2/( E*XI ) ) 

KB2=  CSQRT(  ROH*XB*XH*W** 2/ ( E*XI ) ) 

XK1=  CSQRT( KB 1)*KC*CMPLX( XCA , XSA) 

XK2=-CSQRT(KB1)*KC*CMPLX(XCA,XSA) 

XK3=-CSQRT( KB2)*KC*CMPLX( XCA , XSA ) 

XK4=  CSQRT(KB2)*KC*CMPLX(XCA,XSA) 

KR1=REAL(XK1) 

KR2=REAL(XK2) 

KR3=REAL(XK3) 

KR4=REAL(XK4) 

KI1=AIMAG(XK1) 

KI2=AIMAG(XK2) 

KI3=AIMAG(XK3) 

KI4=AIMAG(XK4) 
C 

C     DEFINED  CONSTANTS  FROM  THE  THREE  REDUCE  PROGRAMS 
C 

P1=V1*XK1 

P2=V1*XK2 

P3=V1*XK3 

P4=V1*XK4 

M 1=EXP( -XL*KI 1 ) * ( COS ( XL*KR 1 ) +V 1*S IN( XL*KR 1 ) ) *XK 1**2 

M2=EXP( -XL*KI 2 ) * ( COS ( XL*KR2 ) +V 1*S IN( XL*KR2 ) ) *XK2**2 

M3=EXP( -XL*KI3 ) *( COS ( XL*KR3 ) +V 1*S IN( XL*KR3 ) ) *XK3**2 


C 

C     IMPEDANCES 
C 


F3=EXP( -XL-KI 3)*( COS ( XL*KR3 ) +V 1*SIN( XL*KR3 ) ) *XK3**: 

F4=EXP(-XL*KI4)*(C0S(XL*I014)+V1*SIN(XL*KR4))*XK4**: 

DEN=(  F4*M3*P2-F4*M3*P1-F4*M2*P3+F4*M2*P1+F4*M1*P3-F4*M1*P2 

-F3*M4*P2+F3*M4*P1+F3*M2*P4-F3*M2*P1-F3*M1*P4+F3*M1*P2 

+F2*M4*P3-F2*M4*P1-F2*M3*P4+F2*M3*?1+F2*M1*P4-F2*M1*P3 

-F1*M4*P3+F1*M4*P2+F 1*M3*P4-F1*M3*P2-F1*M2*P4+F 1*M2*P3 ) 

Al=(  F4*M3*P2-F4*M2*P3-F3*M4*P2+F3*M2*P4+F2*M4*P3-F2*M3*P4) 

/DEN 
A2K-F4*M3*P1+F4*M1*P3+F3*M4*P1-F3*M1*P4-F1*M4*P3+F1*M3*P4) 

/DEN 
A3=(  F4*M2*P1-F4*M1*P2-F2*M4*P1+F2*M1*P4+F1*M4*P2-F1*M2*P4) 

/DEN 
A4=(-F3*M2*P1+F3*M1*P2+F2*M3*P1-F2*M1*P3-F1*M3*P2+F1*M2*P3) 

/DEN 
FF=-V1*E*XI*(A1*XK1**3+A2*XK2**3+A3*XK3**3+A4*XK4**3) 


IMP=2*FF/(V1*W) 

IMPR=REAL(IMP) 
IMPI=AIMAG(IMP) 


SO 


WRITE( 15,1000)  F,IMPR,IMPI 

F=F+20.  0 
100   CONTINUE 

1000   FORMATf IX, D9.  3 ,2X,D12. 5 , 2X,D12. 5) 
CLOSE (UNIT=15) 

STOP 
END 

3.     THE  IMPEDANCE  OF  THE  FINITE  VISCOELASTIC  BEAM  USING 

TIMOSHENKO  BEAM  THEORY 

*  V: 

*  THE  THREE  REDUCE  PROGRAMS  WERE  USED  FOR  EVALUATING  THE  IMPEDANCE  * 
OF  THE  FINITE  VISCOELASTIC  BEAM  USING  TIMOSHENKO  BEAM  THEORY. 

C 

C    THE  FIRST  REDUCE  PROGRAM  DEFINED  SHEAR  FORCE  AND  FOUR  BOUNDARY 

C    CONDITIONS. 

C 

YY:  =Y1*EXP( I*K1*X)+Y2*EXP( I*K2*X)+Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X); 

PP1:  =R1*Y1*EXP( I*K1*X)+R2*Y2*EXP( I*K2*X); 

PP2: =R3*Y3*EXP( I*K3*X)+R4*Y4*EXP( I*K4*X); 

PP:=PP1+PP2; 

MM:=DF(PP,X); 

SS:=-(DF(YY,X)-PP); 

Y0:=SUB(X=0,YY)-Y; 

P0:=SUB(X=0,PP)-P; 

ML:=SUB(X=L,MM); 

SL:=SUB(X=L,SS); 

BYE; 

C 

C    THE  SECOND  REDUCE  PROGRAM  SOLVED  CONSTANT  VARIABLES  (Yl,  Y2 , 

C    Y3  AND  Y4). 

C 

Y0:=Y1+Y2+Y3+Y4-Y; 

P0:=R1*Y1+R2*Y2+R3*Y3+R4*Y4-P; 

ML:  =M1*Y 1+M2*Y2+M3*Y3+M4*Y4; 

SL:=S1*Y1+S2*Y2+S3*Y3+S4*Y4; 

SOLVE(LST(YO,P03ML,SL),Y1,Y2,Y3,Y4); 

BYE; 

C 

C    THE  THIRD  REDUCE  PROGRAM  EVALUATED  THE  SHEAR  FORCE. 

C 

YY:  =Y1*EXP( I*K1*X)+Y2*EXP( I*K2*X)+Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X); 

PP1:  =R1*Y1*EXP( I*K1*X)+R2*Y2*EXP( I*K2*X); 

PP2:  =R3*Y3*EXP(  I*K3*X)+R4*Y4*EXP(  I*K4*X); 

PP:=PP1+PP2; 

SSS:=-KA*GG*AA*(DF(YY,X)-PP); 

SS:=SUB(X=0,SSS); 

BYE; 


*     THIS  FORTRAN  PROGRAM  DEVELOPED  THE  IMPEDANCE  OF  THE  FINITE 
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* 
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SS 

* 
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XI 

* 

XKl, 

* 

XNETA 

•5> 

VI 

VISCOELASTIC  BEAM  USING  TIMOSHENKO  BEAM  THEORY.  * 

DEFINED  VARIABLES; 

COMPLEX  YOUNG'S  MODULUS  OF  THE  VISCOELASTIC  BEAM 

FREQUENCY  (HZ) 

COMPLEX  SHEAR  MODULUS  OF  THE  VISCOELASTIC  BEAM 

GRAVITY  (32$  IN/ SEC 

IMPEDANCE  OF  THE  VISCOELASTIC  BEAM  * 

SHAPE  OF  CROSS  SECTION  * 

FOISSON'S  RATIO  * 

DENSITY  OF  THE  VISCOELASTIC  BEAM 

SHEAR  FORCE  (LBF)  * 

XH,XL  :  DIMENSIONS  OF  THE  BEAM  (WIDTH,  HEIGHT,  LENGTH)      * 
:  AREA  MOMENT  OF  INERTIA  * 

XK2,  XK3,  XK4  :  WAVE  NUMBERS 

:  ENERGY  LOSS  FACTOR  OF  THE  VISCOELASTIC  BEAM 
:  IMAGINARY  (0,1) 

c 

COMPLEX  V1,G,E,XKN1,XKN2,XKD,XK1,XK2,XK3,XK4,R1,R2,R3,R4 

COMPLEX  M1,M2,M3,M4,S1,S2,S3,S4,DEN,A1,A2,A3,A4,SS,IMP 

REAL  F^ARjPOISS^RAVjPIjXB^L^XH.XGjROHjXIjW^XLOG^NETA 

REAL  KR1,KR2,KR3,KR4,KI1,KI2,KI3,KI4,IMPR,IMPI 
C 

OPEN(UNIT=15,FILE='TIM' ) 
C 

C     CONSTANTS 
C 

F=10.  0 

KAR=0.  83 

POISS=0.  45 

GRAV=386.  0 

PI=3.  1415926 

V1=CMPLX(0.  0,1.0) 
C 

C     DIMENSIONS 
C 

XB=1.0 

XL=8. 0 

XH=0. 35 
C 

C     PROPERITIES 
C 

ROH=0. 05505/GRAV 

XI=(XB»XH**3)/12.0 
C 

C     FOUR  DIFFERENT  WAVE  NUMBERS 
C 

DO  100  1=1,100 
W=2-'"PI"F 

XLOG=ABS(LOG10(F/62)) 
XNETA=0. 65*EXP( -0. 52732*XLOG**l. 956) 
C        XNETA=0. 2 
C        XNETA=0.  5 
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XG=0.  00002503*F**3-0.  1752*F**2+457.  5883*F+29280 
C        XG=200000 

G=XG*CMPLX( 1.  0 ,XNETA) 

E=2-G-'-(l+P0ISS) 

XKN1=W**2*( KAR-G+E ) /ROH 

XKN2=  SQRT(  W**4*(  ( KAR*G+E )  /ROH ) **2 
&        -4*KAR*G*E*(W**4-KAR*G*XB*XH*W**2/(R0H*XI)) 
&        /(ROH**2)) 

XKD=2*KAR*  G*E/  ( R0H**2) 

XK1=  CSQRT((XKN1+XKN2)/XKD) 

XK2=-CSQRT((XKN1+XKN2)/XKD) 

XK3=  CSQRT((XKN1-XKN2)/XKD) 

XK4=-CSQRT((XKN1-XKN2)/XKD) 

KR1=REAL(XK1) 

KR2=REAL(XK2) 

KR3=REAL(XK3) 

KR4=REAL(XK4) 

KI1=AIMAG(XK1) 

KI2=AIMAG(XK2) 

KI3=AIMAG(XK3) 

KI4=AIMAG(XK4) 
C 

C     DEFINED  CONSTANTS  FROM  THE  THREE  REDUCE  PROGRAMS 
C 

Rl=( R0H*W**2-KAR*6*XK1**2) / ( V1*KAR*G*XK1) 

R2=(  ROH*W**2-KAR*6*3CK2**2 )  /  (  V1*KAR*G*XK2 ) 

R3=( ROH*W**2 -KAR*G*XK3**2 ) / ( V1*KAR*G*XK3 ) 

R4=( ROH*W**2 -KAR*G*XK4**2 ) / ( V1*KAR*G*XK4 ) 

M1=EXP(-XL*KI1)*(C0S(XL*KR1)+V1*SIN(XL*KR1))*V1*XK1*R1 

M2=EXP ( - XL*KI 2 ) * ( COS ( XL*KR2 ) +V 1 * S IN( XL*KR2 ) ) *V 1 *XK2*R2 

M3=EXP(-XL^KI3)^(COS(XL"KR3)+Vl'vSIN(XL^KR3))-'-Vl-XK3*R3 

M4=EXP( -XL*KI4)*( C0S(XL*KR4)+V1*SIN(XL*KR4) )*V1*XK4*R4 

S 1=EXP( -XL*KT1)*( COS( XL*KR1 )+Vl*SIN( XL*KR1) )*(R1-V1*XK1) 

S2=EXP(-XL-KI2)--(COS(XL^KR2)+Vl"SIN(XL»KR2))-,-(R2-Vl"XK2) 

S3=EXP(-XL*KI3)*(C0S(XL'>KR3)+V1*SIN(XL*KR3))*(R3-V1*XK3) 

S4=EXP(-XL*KI4)*(C0S(XL*KR4)+V1*SIN(XL*KR4))*(R4-V1*XK4) 

DEN=(  S4*M3*R2-S4*M3*R1-S4*M2*R3+S4*M2*R1+S4*M1*R3-S4*M1*R2 
&  - S3*M4*R2+S3*M4*R 1+S3*M2*R4 - S3*M2*R 1 - S  3 *M 1 *R4+S 3*M 1*R2 
&  +S2*M4*R3  -  S  2*M4*R  1  -  S  2  *M3*R4+S  2*M3*R1+S2*M1*R4  -  S2*M  1*R  3 
&        -S1*M4*R3+S 1*M4*R2+S 1*M3*R4-S1*M3*R2 -S1*M2*R4+S 1*M2*R3) 

A 1=(  S4*M3*R2 -S4*M2*R3 - S 3*M4*R2+S3*M2*R4+S 2*M4*R3 - S 2*M3*R4 ) 
&       /DEN 

A2=(-S4*M3*R1+S4*M1*R3+S3*M4*R1-S3*M1*R4-S1*M4*R3+S1*M3*R4) 
&       /DEN 

A3=(  S4*M2*R 1 -S4*M 1*R2 - S 2*M4*R 1+S 2* M 1 * R4+S 1*M4*R2 - S 1*M2*R4 ) 
&       /DEN 

A4=( -S3*M2*R1+S3*M1*R2+S2*M3*R1-S2*M1*R3-S1*M3*R2+S 1*M2*R3 ) 
&       /DEN 

SS=KAR*G*XB*XH*(  A1*(R1-V1*XK1)+A2*(R2-V1*XK2)+A3*(R3-V1*XK3) 
&  +A4*(R4-V1*XK4J  ) 

C 

C      IMPEDANCES 
C 

IMP=2*SS/(V1*W) 

IMPR=REAL(IMP) 

I MP I =A I MAG (IMP) 
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WRITE(15,1000)  F,IMPR,IMPI 

F=F+20. 0 
100   CONTINUE 

1000   FORMATf IX, D9. 3 , 2X,D12.  5 , 2X.D12.  5 ) 
CLOSE (UNIT=15) 

STOP 
END 

4.     THE  IMPEDANCE  OF  THE  FINITE  CONSTRAINED  LAYER  BEAM 

USING  THE  SIXTH  ORDER  BEAM  THEORY 


THE  THREE  REDUCE  PROGRAMS  WERE  USED  FOR  EVALUATING  THE  IMPEDANCE  * 
*    OF  THE  INFINITE  VISCOELASTIC  BEAM  USING  BERNOULLI -EULER  BEAM 
THEORY. 

■)'-  Vc 

c 

C    THE  FIRST  REDUCE  PROGRAM  DEFINED  TRANSVERSE  FORCE  AND  TOO  BOUNDARY 

C    CONDITIONS. 

C 

YY1:  =Y1*EXP( I*K1*X)+Y2*EXP( I*K2*X); 

YY2: =Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X); 

YY3: =Y5*EXP( I*K5*X)+Y6*EXP( I*K6*X); 

UU1: =R1*Y1*EXP( I*K1*X)+R2*Y2*EXP( I*K2*X): 

UU2:  =R3*Y3*EXP(  I*K3*X)+R4*Y4*EXP(  I*K4*X) : 

UU3:  =R5*Y5*EXP( I*K5*X)+R6*Y6*EXP( I*K6*X): 

YYY:=YY1+YY2+YY3; 

UUU: =UU1+UU2+UU3; 

YO:=SUB(X=0,YYY)-Y; 

PO:=SUB(X=0,DF(YYY,X))-P; 

UO:=SUB(X=0,UUU)-U; 

SL:  =SUB(  X=L ,  (  PP*DF( YYY , X ,  3 )+QQ*DF(  YYY ,X)+RR*UUU) ) ; 

ML:=SUB(X=L,(PP*DF(YYY,X,2)+SS*DF(UUU,X))); 

NL:=SUB(X=L,(TT*DF(UUU,X))); 

BYE; 

C 

C    THE  SECOND  REDUCE  PROGRAM  SOLVED  CONSTANT  VARIABLES  (Y4,  Y5 , 

C    AND  Y6). 

C 

YO: =Y1+Y2+Y3+Y4+Y5+Y6-Y; 

P0:=P1*Y1+P2*Y2+P3*Y3+P4*Y4+P5*Y5+P6*Y6-P; 

U0:=R1*Y1+R2*Y2+R3*Y3+R4*Y4+R5*Y5+R6*Y6-U; 

SL:=S1*Y1+S2*Y2+S3*Y3+S4*Y4+S5*Y5+S6*Y6; 

ML:=M1*Y1+M2*Y2+M3*Y3+M4*Y4+M5*YS+M6*Y6; 

NL: =N 1*Y 1+N2*Y2+N3*Y3+N4*Y4+N5*Y5+N6*Y6; 

SOLVE (LST(SL, ML, NL),Y4,Y5,Y6); 

BYE; 

C 

C    THE  THIRD  REDUCE  PROGRAM  SOLVED  VARIABLES  (Yl,  Y2  AND  Y3)  USING 

C    Y4,  Y5  AND  Y6.  THE  THIRD  REDUCE  PROGRAM  EVALUATED  THE  SHEAR  FORCE. 

C 

Y4 

Y5 

Y6 

YO 


=B 1*Y 1+C 1*Y2+D 1*Y3; 
=B2*Y1+C2*Y2+D2*Y3; 
=B3*Y1+C3*Y2+D3*Y3; 

=Y1+Y2+Y3+Y4+Y5+Y6-Y; 
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P0:=P1*Y1+P2*Y2+P3*Y3+F4*Y4+P5*Y5+P6*Y6-P; 

U0:=R1*Y1+R2*Y2+R3*Y3+R4*Y4+R5*Y5+R6*Y6-U; 

S0LVE(LST(Y0,P0,U0),Y1,Y2,Y3); 

YYY1:  =Y1*EXP( I*K1*X)+Y2*EXP( I*K2*X)i 

YYY2: =Y3*EXP( I*K3*X)+Y4*EXP( I*K4*X) 

YYY3:  =Y5*EXP( I*K5*X)+Y6*EXP( I*K6*X): 

YYYY: =YYY1+YYY2+YYY3; 

UUU1: =R1*Y1*EXP( I*K1*X)+R2*Y2*EXP( I*K2*X); 

UUU2: =R3*Y3*EXP( I*K3*X)+R4*Y4*EXP( I*K4*X); 

UUU3:  =R5*Y5*EXP( I*K5*X)+R6*Y6*EXP( I*K6*X) ; 

UUUU: =UUU1+UUU2+UUU3; 

S  S  S :  =P  P*DF  (  YYYY ,  X ,  3 )  +QQ*DF  ( YYYY ,  X )  +RR*UUUU ; 

SS:=SUB(X=0,SSS); 

BYE; 


* 

* 
* 

* 

-V 

* 
* 
* 

■k 
* 
* 

* 

* 


THIS  FORTRAN  PROGRAM  DEVELOPED  THE  IMPEDANCE  OF  THE  CONSTRAINED  * 
LAYER  BEAM  USING  THE  6TH  ORDER  BEAM  THEORY.  * 


DEFINED  VARIABLES; 

El,  E3  :  COMPLEX  YOUNG'S  MODULUS  OF  THE  CONSTRAINED  LAYER 

F  :  FREQUENCY  (HZ) 

G  :  COMPLEX  SHEAR  MODULUS  OF  THE  VISCOELASTIC 

GRAV  :  GRAVITY  (386  IN/SEC**2) 

HI,  H3  :  HEIGHT  OF  THE  CONSTRAINED  LAYER 

H2  :  HEIGHT  OF  THE  VISCOELASTIC 

II,  12  :  AREA  MOMENT  OF  INERTIA 

IMP  :  IMPEDANCE  OF  THE  CONSTRAINED  LAYER  BEAM 

KAR  :  SHAPE  OF  CROSS  SECTION 

POISS  :  POISSON'S  RATIO 

ROH1  :  DENSITY  OF  THE  CONSTRAINED  LAYER 

ROH2  :  DENSITY  OF  THE  VISCOELASTIC 

SS  :  SHEAR  FORCE  (LBF) 

XB,  XL  :  DIMENSIONS  OF  THE  BEAM  (WIDTH,  LENGTH) 

XK1,  XK2,  XK3,  XK4,  XK5 ,  XK6  :  WAVE  NUMBERS 

XNETA  :  ENERGY  LOSS  FACTOR  OF  THE  VISCOELASTIC  BEAM 

VI  :  IMAGINARY  (0,1) 


* 


if 


-y?VcycyryoV-.vyoV}Viv-Wr-w?ycyc-Vivyw-ycVc-v^ 


COMPLEX*32  VI , G , BB , PI , P2 , P3 , P4 , P5 , P6 ,R 
COMPLEX*32  S 1 , S2 , S3 , S4 , S5 , S6 ,M1 ,M2 , M3 , 
COMPLEX*32  B1,C1,D1,B2,C2,D2,B3,C3,D3, 
COMPLEX  B,D,XK1,XK2,XK3,XK4,XK5,XK6 

REALM6  F,PI,GRAV,XB,XL,XG,H1,H2,H3,E1 
REAL*16  CC,XM,T1(30),T2(30),T3(30),Q1( 
REAL*16  XL0G,W,Z,C,BR,BI,CR,DR,DI,FT1, 
REAL''- 16  FT2,FQ2,DTT2,DTQ2,DQQ2,DQT2,FT 
REAL* 16  XKR1,XKR2,XKR3,XKI1,XKI2,XKI3, 
REAL--'"  16  KI1,KI2,KI3,KI4,KI5,KI6,PP,S3, 
REAL  XNETA 

OPEN(UNIT=15,FILE='SAN' ) 


1,R2,R3,R4,R5,R6,QQ,RR 

M4,M5,M6,N1,N2,N3,N4,N5,N6 

DEN1,DEN2,A1,A2,A3 

,XXK1,XXK2,XXK3,IMP 

,E3,I1,I3,DD,R0H1,R0H2 

30),Q2(30),Q3(30) 

FQ1 , DTT1 , DTQ 1 . DQQ 1 , DQT1 

3 , FQ3 , DTT3 , DTQ3 , DQQ3 , D0T3 

KR1 , KR2 , KR3 , KR4 , KR5 , KR6 

TT,IMPR,IMPI 


S5 


c 

G 


C 

c 
c 

c 


CONSTANTS 

F=2010 

PI=3. 1415926 
GRAV=386.  0 
V1=CMPLX(0.  0,1.  0) 

DIMENSIONS 

XD=1.  0 
XL=8. 0 
Hl=0.0625 
H2=0. 065 
H3=0.  0625 

PROPERITIES 

El=10000000.  0 

E3=10000000.  0 

I1=XB*H1**3/12. 

I3=XB*H3**3/12. 

DD=E1*I1+E3*I3 

R0H1=0.  09832/GRAV 

R0H2=0.  03663/GRAV 

CC=Hl/2. +H2+H3/2. 

XM=XB*( ( H1+H3 ) *R0H1+H2*R0H2 ) 

INITIAL  GUESS  WAVE  NUMBERS 

Tl(l)=-4.  31911862 
T2(l)  =  -0.  05139600 
T3(l)=  3.  77697017 
Ql(l)  =  -0.  33835599 
Q2(l)=-0.  05104052 
Q3(l)=-0.  27669665 

FIND  SIX  DIFFERENT  WAVE  NUMBERS  USING  NEWTON'S  METHOD 
SIX  DIFFERENT  WAVE  NUMBERS  ARE  XIU ,XK2 ,XK3,XK4,XK5 ,XK6 

DO  100  1=1,100 
W=2*PI*F 

XLOG=ABS(LOG10(F/800)) 
XNETA=1.  1*EXP(  -0.  3906*XLOG**1.  53) 
XNETA=0. 5 
XNETA=1.  0 

XG=-0.  00001258*F**2+0.  2472*F+74.  988 
XG=1000 

G=XG*CMPLX( 1.  0,XNETAj 
BB=G*( E 1*H1+E3*H3 ) / ( H2*E 1*H1*E3*H3 ) 
Z=cC--» 2*E  1  "H1--E3-H3/ ( DD*( E i*Hl+E3*H3 ) ) 
B=BB*(1.  O+Z) 
C=-XM*W**2/DD 
D=-BB*XM*W**2/DD 
BR=REAL(B) 
BI=AIMAG(B) 
CR=C 
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DR=REAL(D) 

DI=AIMAG(D) 
DO  200  J=l,15 

FTl=Tl(J)**3+BR*Tl(J)**2+(-3*Ql(J)**,-2-2*BI'VQl(J)+CR)',fTl(J) 
&  +(-BR*Ql(J)**2+DR) 

FQl=Ql(J)**3+BI*Ql(J)**2+(-3*Tl(J)**2-2*BR*Tl(J)-CR)*Ql(J) 
&  +(-BI*Tl(J)**2-DI) 

DTT1=3*T1( J)**2+2*BR*T1( J) 
&  +(  - 3*Q 1  ( J ) ** 2  - 2*B I*Q1(  J)+CR) 

DTQ1K-6*T1(J)-2*BR)*Q1(J)-2*BI*T1(J) 
DQQ 1=3*Q  1  ( J )  **2+2*BI  *  Q 1  ( J ) 
&  +(-3*Tl(J)**2-2*BR*Tl(J)-CR) 

DQT1<-6*Q1(J)-2*BI)*T1(J)-2*BR*Q1(J) 

Tl( J+l )=T1( J) -(  DQQ1*FT1-DTQ1*FQ1 )/(DTTl*DQQl-DTQl*DQTl) 
Ql( J+l )=Q1( J) - ( -DQT1*FT1+DTT1*FQ1) / ( DTT1*DQQ1-DT01*D0T1 ) 
FT2^T2(J)**3+BR*T2(J)**2+(-3*Q2(J)**2-2*BI*Q2(J)+CR)*T2(J) 
&  +(-BR*Q2(J)**2+DR) 

FQ2^2(J)^3+BI*Q2(J)**2+(-3*T2(J)**2-2*BR*T2(J)-CR)*Q2(J) 
5c  +(-BI*T2(J)**2-DI) 

DTT2=3*T2  ( J )  **2+2*BR*T2  ( J ) 
5c  +(-3*Q2(J)**2r.2*BI*Q2(J)+CR) 

DTQ2K-6*T2(J)-2*BR)*Q2(J)-2*BI*T2(J) 
DQQ2=3*Q2(J)**2+2*BI*Q2(J) 
5c  +(  -3*T2(  J)**2-2*BR*T2(  J)  -CR) 

DQT2=(-6*Q2(J)-2*BI)*T2(J)-2*BR*Q2(J) 

T2(J+l)=T2(J)-(  DQQ2*FT2-DTQ2*FQ2)/(DTT2*DQQ2-DTQ2*D0T2) 
Q2(J+l)=Q2(J)-(-DQT2*FT2+DTT2*FQ2)/CDTT2'>DQQ2-DT02',-DQT2) 
FT3=T3(J)**3+BR*T3( J)**2+( -3*Q3( J)**2-2*BI*Q3( J)+CR)*T3( J) 
5c  +(-BR*Q3(J)**2+DR) 

FQ3=Q3(J)**3+BI*Q3(J)**2+(-3*T3(J)**2-2*BR*T3(J)-CR)*Q3(J) 
&  +(-BI*T3(J)**2-DI) 

DTT3=3*T3 ( J ) **2+2*BR *T3 ( J ) 
5c  +(  -3*Q3(  J)**2-2*BI*Q3(  J)+CR) 

DTQ3=(-6*T3(J)-2*BR)*Q3(J)-2*BI*T3(J) 
DQQ3=3*Q3(J)**2+2*BI*Q3(J) 
&  +( -3*T3( J)**2-2*BR*T3( J) -CR) 

DQT3=(-6*Q3(J)-2*BI)*T3(J)-2*BR*Q3(J) 

T3( J+l )=T3( J) -(  DQQ3-FT3 -DTQ3*FQ3 ) / ( DTT3*DQQ3 -DTQ3*DQT3 ) 
Q3( J+l )=Q3( J) -( -DQT3*FT3+DTT3*FQ3 ) / ( DTT3*DQQ3 -DTQ3*DQT3 ) 
200    CONTINUE 

XKR1=  T1(J) 

XKR2=  T2(J) 

XKR3=  T3(J) 

XKI1=  Q1(J) 

XKI2=  Q2(J) 

XKI3=  Q3CJ) 

XXK 1 =C  MP  LX ( XKR 1 , XK 1 1 ) 

XXK2=CMPLX( XKR2 , XKI2 ) 

XXK3=CMPLX( XKR3 , XKI 3 ) 

XK1=  CSQRT(XXKl) 

XK2=-CSQRT(XXK1) 

XK3=  CSQRT(XXK2) 

XK4=-CSQRT(XXK2) 

XK5=  CSQRT(XXK3) 

XK6=-CSQRT(XXK3) 

KR1=  REAL(XKl) 
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KR2=  REALCXK2) 

KR3=  REAL(XK3) 

KR4=  REAL(XK4) 

KR5=  REAL(XK5) 

KR6=  REAL(XK6) 

KT1=AIMAG(XK1) 

KI2=AIMAG(XK2) 

KI3=AIMAG(XK3) 

KI4=AIMAG(XK4) 

KI5=AIMAG(XK5) 

KI6=AIMAG(XK6) 
C 

C     DEFINED  CONSTANTS  FROM  THE  THREE  REDUCE  PROGRAMS 
C 

P1=V1*XK1 

P2=V1*XK2 

P3=V1*XK3 

P4=V1*XK4 

P5=V1*XK5 

P6=V1*XK6 

Rl=Vl*BB*Z*DD*XKl/( ( XK1**2+BB)*E3*H3*CC) 

R2=V1*BB*Z *DD*XK2/ ( ( XK2**2+BB )*E3*H3*CC ) 

R3=Vl-vBB^Z"DD»XK3/((XK3-'""-2+BB)''>-E3^H3^CC) 

R4=V  1*BB*Z*DD*XK4/  ( ( XK4**2+BB )  *E3*H3 *CC ) 

R5=Vl*BB*Z*DD*XK5/( ( XK5**2+BB)*E3*H3*CC ) 

R6=Vl*BB*Z*DD*XK6/( f XK6**2+BB)*E3*H3*CC ) 

PP=DD 

QQ=-DD*BB*Z 

RR*DD*BB*Z*(E1*H1+E3*H3 ) / ( CC*E 1*H1) 

SS=-CC*E3*H3 

TT=E3*H3 

S 1=( - V1*PP*XK1**3+V 1 *QQ*XK 1+RR*R 1 ) -EXP ( - XL *K I 1 ) 
&   *(C0S(XL*KR1)+V1*SIN(XL*KR1)) 

S2=( -V1*PP*XK2**3+V1*QQ*XK2+RR*R2)*EXP( -XL*KI2) 
&   *( COS ( XL*KR2 ) +V 1*S IN( XL*KR2 ) ) 

S3=(-V1*PP*XK3**3+V1*QQ*XK3+RR*R3)*EXP(-XL*KI3) 
&   *(COS(XL*KR3)+Vl*SIN(XL*KR3)) 

S4=(-V1*PP"XK4**3+V1*QQ*XK4+RR*R4)*EXP(-XL"KI4) 
&   *(C0S(XL*KR4)+V1*SIN(XL*KR4) ) 

S5=( -V1*PP*XK5**3+V1*QQ*XK5+RR*R5)*EXP( -XL*KI5 ) 
&   *(  COS  ( XL*KR5 )  +V  1*SIN(  XL*KR5 ) ) 

S6=(  -V1*PP*XK6**3+V1*QQ*XK6+RR*R6 )*EXP(  -XL--KI6 ) 
&   *(C0S(XL*KR6)+V1*SIN(XL*KR6) ) 

Ml=( V1*SS*R1*XK1-PP*XK1**2)*EXP( -XL*KI1 ) 
&   * ( C  0  S ( XL*KR 1 ) +V 1 *  S I N ( XL*KR 1 ) ) 

M2=( V 1*SS*R2*XK2 -PP*XK2**2 ) -EXP ( - XL- K I 2 ) 
&   *(  COS  ( XL*KR2 )  +V 1*S IN(  XL*KR2) ) 

M3KV1*SS*R3*XK3-PP*XK3**2)*EXP(-XL*KI3) 
&   *(C0S(XL*KR3)+V1*SIN(XL*KR3)} 

M4=(  V1*SS*R4*XK4-PP*XK4**2 )*EXP(  -XL-KI4) 
&   *( COS ( XL*KR4 ) +V 1*SIN( XL*KR4 ) ) 

M5KV1*SS*R5*XK5-PP*XK5**2)*EXPC-XL*KI5) 
&   *(  COS ( XL*KR5 )  +V  1*S IN(  XL*KR5 ) ) 

M6=( V1*SS*R6*XK6 -PP*XK6**2) *EXP( -XL*KI6 ) 
Sc        *(  COS ( XL*KR6 )  +V  1*SIN(  XL*KR6 )  ) 

N 1=V1*TT*R 1*XK1*EXP(  -XL*KI 1 ) 
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&       *(C0S(XL*KR1)+V1*SIN(XL*KR1)) 

N2=V1*TT*R2*XK2*EXP( -XL*KI2) 
&       *(C0S(XL*KR2)+V1*SIN(XL*KR2)) 

N  3= V  1*TT*  R3*XK3  *  EXP  (  -  XL*  K 1 3 ) 
&       *('  COS ( XL*KR3 )  +V1*S IN(  XL*KR3 ) ) 

N4=V i*TT*R4*XK4*EXP ( - XL*K I 4 ) 
&        *( COS( XL*KR4)+V1*SIN( XL*KR4) ) 

N5=V1*TT*R5*XK5*EXP( -XL*KI5 1 
&        *( COS( XL*KR5 )+Vl*SIN( XL*KR5 ) ) 

N6=V1*TT*R6*XK6*EXP( -XL*KI 6 ) 
&        *(C0S(XL*KR6)+V1*SIN(XL*KR6) ) 

B 1=( -H6*M5*S 1+N6*M1*S5+N5*M6*S 1 -N5*M1*S6 -N1*M6*S5+N1*M5*S6 ) 
6c        /(N6*M5*S4-N6*M4*S5-N5*M6*S4+N5*M4*S6+N4*H6*S5-N4*M5*S6) 

Cl=(-N6*M5*S2+N6*M2*S5fN5*M6*S2-N5*M2*S6-N2*M6*S5+N2*iM5*S6) 
*        /  ( N6*M5*S4-N6*M4*S5  -N5*M6*S4+N5*M4*S6+N4*M6*S5  -N4*U5*S6 ) 


& 

C2=(N6*M4*S2-N6*M2*S4-N4*M6*S2+N4*M2=>S6+N2*M6*S4-N2*M4*S6) 
&        /(N6*M5*S4-N6*M4*S5-N5*M6*S4+N5*M4*S6+N4*tl6*S5-N4*M5*S6) 

D2=( N6*M4*S3 -N6*M3*S4 -N4*M6*S3+N4*M3*S6+N3*M6*S4-N3*M4*S6 ) 
&        /(N6*M5*S4-N6*M4*S5-N5*M6*S4+N5*M4*S6+N4*M6*S5-N4*M5*S6) 

B3=( -N5*M4*S 1+N5*M 1*S4+N4*M5*S 1 -N4*M 1*S5 -N 1*M5*S4+N 1*M4*S5 ) 
&        /(N6*M5*S4-N6*M4*S5-N5*M6*S4+N5*M4*S6+N4*M6*S5-N4*M5*S6) 

C3=( -N5*M4*S2+N5*M2*S4+N4*M5*S2 -N4*M2*S5 -N2*M5*S4+N2*M4*S5 ) 
6c 

D3= 

DEN1=( R6*P5 -R5*P6 )*( D3*C2*B 1+D3*C2 -D3*B2*C 1 -D3*B2 -C3*D2*B 1 
6c  -C3*D2+C3*B2*D1+C3*B2+B3*D2*C  1+B3*D2 

5c  -B3*C2*D1-B3*C2) 

5c  +(  R4*P6  -R6*P4 )  *(  D3*C2*B  1  -D3*B  2*C  1  -D3*C  1+D3*B  1  -C3*D2*B  1 
5c  +C3*B2*D1+C3*D  1  -C3*B  1+B3*D2*C  1  -B3*C2*D1 

6c  -B3*D1+B3*C1) 

6c  +( R5*P4-R4*P5 )  *(  D3*C2*B  1  -D3*B2*C  1 -C3*D2*B  1+C3*B2*D  1 
6c  +B3*D2*C1-B3*C2*D1+D2*C1-D2*B1-C2*D1 

6c  +C2*B1+B2*D1-B2*C1) 

DEN2=(R3*P6-R6*P3)*(C3*B2+G3*B1+G3-B3*C2-B3*C1-B3) 

&  +(R6*P2-R2*P6)*(D3*B2+D3*B1+D3-B3*D2-B3*D1-B3) 

6c  +(R1*P6-R6*P1)*(D3*C2+D3*C1+D3-C3*D2-C3*D1-C3) 

&  +(R5*P3-R3*P5)*(C3*B2-B3*C2-C2*B1-C2+B2*C1+B2) 

6c  +(R2*P5-R5*P2)*(D3*B2-B3*D2-D2*B1-D2+B2*D1+B2) 

6c  +(  R5*P  1  -R1*P5  )*(  D3*C2  -C3*D2  -D2*C  1  -D2+C2*D1+C2  ) 

&  +( R4*P3 -R3*P4)*( C3*B 1 -B3*C 1+C2*B 1 -B2*C 1 -C 1+B 1 ) 

&  +(  R2*P4  -R4*P2  )  *(  D3*B  1  -  B  3*D  1+D2*B  1  -  B  2*D  1  -D  1+B  1 ) 

6c  +( R4*P  1  -R1*P4  )*(  D3*C  1  -C3*D  1+D2*C  1  -C2*D  1  -D 1+C 1 ) 

6c  +(R3*P2-R2*P3)*(B3+B2+B1+1) 

6c  +(R1*P3-R3*P1)*(C3+C2+C1+1) 

&  +(R2*P1-R1*P2)*(D3+D2+D1+1) 

Al=(    R6*(    P5*D3*C2-P5*C3*D2+P4*D3*C1-P4*C3*D1-P3*G3+P2*D3) 

5c  +  R5*( -P6*D3*C2+P6*C3*D2+P4*D2*C1-P4*C2*D1-P3*C2+P2*D2) 

6c  +  R4*(  -P6*D3*C1+P6*C3*D1-P5*D2*C1+P5*C2*D1-P3*C1+P2*D1) 

oc  +  R3*(    P6*C3+P5*C2+P4*C1+P2) 

5c  +  R2*(-P6*D3-P5*D2-P4*D1-P3)    )/(DENl+DEN2) 

A2=(    R6*(-P5*D3*B2+P5*B3*D2-P4*D3*B1+P4*B3*D1+P3*B3-P1*D3) 


S9 


&     +  R5*(  P6*D3*B2-P6*B3*D2-P4*D2*B1+P4*B2*D1+P3*B2-P1*D2) 
&    +  R4*(  P6*D3^-P6*B3*D1+P5*D2*B1-P5*B2*D1+P3*B1-P1*D1) 
6c    +  R3*(-P6*B3-P5*B2-P4*B1-P1) 
6c    +  Rl*(  P6*D3+P5*D2+P4*D1+P3)  )/(DENl+DEN2) 

A3=(  R6*(  P5*C3*B2-P5*B3*C2+P4*C3*B1-P4*B3*C1-P2*B3+P1*C3) 
&    +  R5*(  -P6*C3*B2+P6*B3*C2+P4*C2*B1-P4*B2*C1-P2*B2+P1*C2) 
6c    +  R4*(  -P6*C3*B1+P6*B3*C1-P5*C2*B1+P5*B2*C1-P2*B1+P1*C1) 
6c    +  R2*(  P6*B3+P5*B2+P4*B1+P1) 
6c    +  R1*(-P6*C3-P5*C2-P4*C1-P2)  )/(DENl+DEN2) 

C 

C      IMPEDANCES 

G 

IMP=2*(  ( R4*RR+V1*QQ*XK4-V1*PP*XK4**3 ) *(D1*A3+C  1*A2+B  1*A1 ) 
&    +   ( R5  *RR+V 1*QQ*XK5 - V 1*PP*XK5 ** 3 ) * ( D  2* A3+G 2* A2+B 2* A 1 ) 
6c    +   ( R6*RR+V1*QQ*XK6  -V1*PP*XK6**3 )  *(  D3*A3+C3*A2+B3*A1 ) 
&    +  RR*(R3*A3+R2*A2+R1*A1) 

&    +  V1*QQ*(XK3*A3+XK2*A2+XK1*A1) 

&    -V1*PP*(  K3**3*A3+K2**3*A2+K1**3*A1)  )  /  ( V1*W) 

IMPR=REAL(IMP) 

IMPI=AIMAG(IMP) 

WRITE( 15 , 1000)F, IMPR, IMPI 

F=F-20 

T1(1)=T1(J) 

T2(1)=T2(J) 

T3(1)=T3(J) 

Q1(1)=Q1(J) 

Q2(1)=Q2(J) 

03(1)=Q3(J) 
100    CONTINUE 
1000   FORMAT( IX, D9. 3 , 2X,D12. 5 , 2X,D12. 5) 

CL0SE(UNIT=15) 

STOP 

END 
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APPENDIX  B. 
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Figure  49.       Real  part  of  the  20"  viscoelastic  beam  waveguide  absorber  impedances 
at  the  center  of  the  beam. 
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Figure  50.       Imaginary  part  of  the  2U"  viscoelastic  beam  waveguide  absorber  impe- 
dances at  the  center  of  the  beam. 
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Figure  51.  The  driving  point  frequency  response  of  the  test  plate  nith  a  16"  viscoe- 
lastic  heam  waveguide  absorber  (dashed)  at  location  1  and  without 
(solid). 
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Figure  52.      Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  10"  viscoe- 
lastic  beam  waveguide  absorber  at  location  I  and  without. 
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Figure  53.      The  driving  point  frequency  response  of  the  test  plate  with  a  20"    con- 
strained  layer   beam   waveguide  absorber  (dashed)   at  location    1    and 

without  (solid). 
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Figure  54.       Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  20"  con- 
strained laver  beam  waveguide  absorber  at  location  1  and  without. 
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Figure  55.  The  driving  point  frequency  response  of  the  test  plate  with  a  16"  viscoe- 
Iastic  beam  waveguide  absorber  (dashed)  at  location  2  and  without 
(solid). 


97 


2 

W 

or 

ft: 

09 

> 

g 

< 

O 


?li.fli. 


o 

IMS   • 


a 

-a. 


-•.a; 


•  - 


a  : 


"<• 0" 


I  I  I    I    I     I       I 


X 

Ed 

03 
XX 
WO 

a« 
o«< 

§§* 

all 

ii  ii 

a  • 


T-rr 


B 


z 

Ed 

Or 
U 

OS 
b. 


Houvi  3Nidwva  ivaow 


4-o 


JOW 


Figure  56.      Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  16"  viscoe- 
lastie  beam  waveguide  absorber  at  location  2  and  without. 


Figure  57.  The  driving  point  frequency  response  of  the  test  plate  with  a  16"  con- 
strained layer  beam  waveguide  absorber  (dashed)  at  location  2  and 
without  (solid). 
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Figure  58.       Modal  damping  factors  vs.  frequency  of  Hie  test  plate  with  a  16''  con- 
strained laver  beam  waveguide  absorber  at  location  2  and  without. 
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Figure  59.  The  driving  point  frequency  response  of  the  test  plate  with  a  20"  con- 
strained layer  beam  waveguide  absorber  (dashed)  at  location  2  and 
without  (solid). 
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Figure  60.       Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  20"  con- 
strained layer  beam  waveguide  absorber  at  location  2  and  without. 
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Figure  61.  The  driving  point  frequency  response  of  the  test  plate  with  a  16"  viscoe- 
lastic  beam  waveguide  absorber  (dashed)  at  location  3  and  without 
(solid). 
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Figure  62.      Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  16"  viscoe- 
lastic  beam  waveguide  absorber  at  location  3  and  without. 


104 


Figure  63.  The  driving  point  frequency  response  of  the  test  plate  with  a  16"  con- 
strained layer  beam  waveguide  absorber  (dashed)  at  location  3  and 
without  (solid). 
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Figure  64.       Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  16"  con- 
strained laver  beam  waveguide  absorber  at  location  3  and  without. 
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Figure  65.  The  driving  point  frequency  response  of  the  test  plate  with  a  20"  con- 
strained layer  beam  waveguide  absorber  (dashed)  at  location  3  and 
without  (solid). 
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Figure  66.      Modal  damping  factors  vs.  frequency  of  the  test  plate  with  a  20"  con- 
strained layer  beam  waveguide  absorber  at  location  3  and  without. 
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